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(54) Turbo-reception method and turbo-receiver for a MIMO system 



(57) An impluse response h mn (q) of each transmis- 
sion path is estimated from N received signals r m (m=1 , 
M) and a known signal (for a number of users equal 
to N, n=1, ■•. N). MxN matrix H (q) having h mn (q) as an 
element and a QxQ matrix H having H(q) as an element 
are determined (where Q represents a number of mul- 
tipaths of each transmitted wave and q=0, Q-1). A 
soft decision value b' n (k) is determined from decoded 
^2 [b n (k)], and this is used to generate an interference 
component matrix B'(k) to generate an interference rep- 
lica H*B*(k). The interference replica H-B'(k) is subtract- 
ed from a received matrix y(k) to determine y*(k). y(k) 
and H are used to determine an adaptive filter coefficient 
w n (k) to be applied to an n-th user in order to eliminate 
residual interference components in y*(k) according to 
the minimum mean square error criteria. y(k) is passed 
through w n (k) to provide a log-likelihood ratio as a re- 
ceived signal from the user n from which interferences 
are eliminated. 
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D s riptl n 

BACKGROUND OF THE INVENTION 

[0001] The invention relates to turbo-reception method and a turbo-receiver as may be used in a mobile communi- 
cation, for example, and which apply an iterative equalization utilizing a turbo-coding technique to waveform distortions 
which result from interferences. 

[0002] A task in the mobile station communication business is how to construct a system capable of acquiring a 
multitude of users on a limited frequency domain with a high quality. A multi-input multi-output (MIMO) system is known 
in the art as means for solving such a task. The architecture of this system is shown in Fig. 30A where a plurality of 
transmitters S1 to SN transmit symbols Cl (i) to c N (i) at the same time and on the same frequency, and the transmitted 
signals are received by an MIMO receiver equipped with a plurality of antennas #1 to #M. The received signals are 
processed by the receiver, which esjjmates A transmitted symbols d(!) to ctfl) from the respective transmitters S1 to 
SN and delivers them separately as c^i) to c N (i) to output terminals Out 1 to Out N. 

[0003] Up to the present time, a study of a specific implementation of an MIMO receiver in an MIMO system is not 
yet satisfactonly warranted. If one attempts to construct an MIMO receiver in an MIMO system on the basis of MLSE 
(maximum likelihood estimation) criteria, denoting the number of transmitters by N and the number of multi-paths 
through which a wave transmitted from each transmitter reaches the MIMO receiver by Q, the quantity of calculation 
required for the MIMO receiver will be on the order of 2<Q-DN and will increase even more voluminously with an increase 
in the number of transmitters N and the number of multi-paths Q. When information from a single user is transmitted 
as parallel signals, which are then received, a separation of individual parallel signals from each other requires a 
quantity of calculation, which increases exponentially with number of multi-paths. Accordingly, the present invention 
proposes herein an improved calculation efficiency turbo-reception method for a plurality of channel signals To start 
with, an existing turbo-receiver for a single user (single transmitter) or a single channel transmitted signal, which illus- 
trates the need for the present invention, will be described. 

Turbo-receiver for single user 

[0004] An exemplary arrangement for a transmitter and a receiver is illustrated in Fig. 31 . In a transmitter 1 0 infor- 
mation senes c(i) is encoded in encoder 1 1 , and an encoded output is interleaved (or rearranged) by an interlea'ver 1 2 
before it is input to a modulator 1 3 where it modulates a carrier signal, the resulting modulated output being transmitted 
The transmitted signal is received by a receiver 20 through a transmission path (each channel of multipath) In the 
receiver 20, a soft input soft output (SISO: single-input single output) equalizer 21 performs an equalization of delayed 
waves. At the input to the equalizer 21, the received signal is generally converted into a baseband, and the received 
baseband signal is sampled with a frequency which is equal to or greater than the frequency of symbol signals of 
information senes in the transmitted signal to be converted into a digital signal, which is then input to the equalizer 21 
[0005] For the single user, this corresponds to N=1 in Fig. 30A. and a received output form each reception antenna 
ffm in (m=1 , 2, • •-, M) can be represented as follows: 

rmW =2 q=0 Q ' 1 h m (q) b(k-q) + v m (k) (1) 

where m represents an antenna index, h a channel value (a transmission path impulse response: a transmission path 
charactenstic), b(k-q) a transmitted symbol from a user (transmitter 1), and v m (k) an internal thermal noise of the 
receiver 20. All outputs from the antennas #1 to #M are denoted by a matrix as indicated by an equation (2) to define 
an equation (3). ^ w 

r(k)=[r 1 (k)r 2 (k)...r M (k)] T (2) 
=Z<fO°'' H(q) b(kHq) + v(k) (3) 

where 



v(k) =[v 1 (k)v 2 (k)...v M (k)] T 



(4) 
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10 



H(q)=[h 1 (q)...h M (q)] T 



(5) 



Itistob not d that Q T represents an inverted matrix. In cosideration of the number of channels Q of the mutipath the 
5 following matrixes and matrix are defined: 



where 



y(k)^[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 



=Hb(k)+n(k) 



(6) 
(7) 



15 



20 



H = 



Ho 



0) 



0 



H 



(0-1) 



'(0) 



0 



H 



(0-1). 



(8) 



25 



b(k-q)=[b(k+Q-1)b(k+Q-2)...b(k-Q+1)] T 



(9) 



30 



n(k) =[v T (k+Q-1 )v T (k+Q-2)...v T (k)] T 



(10) 



[0006] r(k) as defined above is input to the SISO equalizer 21. which is a linear equalizer, deriving a log-likelihood 
ratjo A,(LLR) of a probability that each encoded bit {b(i» is equal to +1 to a probability that it is -1 as an equalization 
output. 



35 



A 1 [b(k)]=log 



Prtb(k) = -1|y(k)] 



(11) 



40 



45 



=Xi[b(k)]+*P[b(k)] 



(12) 



where X,[b(k)] represents an extrinsic information fed to a succeeding decoder 24 and X p [b(k)] a priori information 
applied to the equalizer 21 . The log-likelihood ratio A^bfk)] is fed to a subtracter 22 wherl the a priori information X, 
[b(k)] is subtracted therefrom. The result is then fed through a deinterleaver 23 to an SISO channel decoder 24 which 
calculates a log-likelihood ratio A 2 as follows: 



A 2 [b(i)] = log 



Prtb(i) = -t-1|X 1 [b(i)] i i = 0,- l B-1] 
Pr[b(i) = -1|X 1 [b(i)],i = 0, -,B-1] 



50 



B : frame length 



(13) 



- A*[b(i)]+X?[b(i)] (14) 

55 ^ Te ^F)}. represents an e^nsic information which is applied as ^rb(k)] to the equalizer 21 during the iteration 
while Mbfk)] is applied as a priori information \ F [b(i)] to the decoder 24. In a subtracter 25, X,M)] is subtracted from 
A 2 [b(,)]. and the result is fed through an interieaver 26 to the equalizer 21 and the subtracter 22. In this manner the 
qualization and the decoding are iterat d to achieve an improvement of an rror rat . 
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[0007] Todescrib th prestage qualizer21 in detail, the calculation of a linear filt rreso nseannn ritoaro^n 
d6SC " b ^singtheaprioninfonnation^pXkHforth ^r" ^ZK^SSSL: 

b'(k)=tan h[x£ [b(k)]/2] (15) 

is calculated. Using the estimate and a channel matrix H, an interference component or a replica H-b'(k) of the inter- 
ference component is reproduced and subtracted from the received signal. Thus. 

y'(k)=y(k)-Hb'(k) (16) 

=H(b(k)-b'(k))+n(k) (17) 

where, 

b-(k)=[b , (k+Q-1)...0...b'(k-Q+1)] T (18) 

Because the replica H-b'(k) of the interference component cannot be always a correct replica, the interference com- 
residl^L nt e rf COmP ' ete,y e ' iminated bV eqUati ° n (16) - So 8 ,ine3r f ' ,l3r coefficient P w(k SSSSL^i 

EZ1 fSZ component determined aocordin9 10 the MMSE (minimum mean i-» *™> 



w(k)=arg min ||w H (k)y(k)-b(k)|| 2 (19) 

30 ££" De2mo^knntfl?„ n8P08 ! i0n T! ! |a : 0ml - WM WhiCh mimimi2eS 013 3 < U3tion < 19 > is determined. 
E2L£! £n k ^ T 6 ' S descnbed ln Daf y and Xiandong Wang. Tow Complexity Turbo- 

SnffS h t 6rS ^ Channe ' S (^'eetarnu.edu/Reynolds/). A major achievement of this technique lies l a 
significant reduction m the quantity of calculation. The quantity of calculation according to a conventional MLSE tJroo 
has been proportional to the order of *" while a suppression to the order of Q3 , 8 enabled by this Snique t wiN 

35 ? See h " h a H^ ( t k) / (k) rSpreSentS a " ° UtpUt fmm "» 3 « uali2er 21 • a " d is to calculate^ ,[b(k) 'which is then fed 
hrough the de.nterteaver 23 to the decoder 24 to be used in the decoding calculation. 

EL PUrP ° 8e °! , equa,ization in ^ uali2er 21. it is necessary to estimate the channel value (transmission 
path impulse response) h appearing in the equation (1 ). This estimation is hereafter referred to as a channel esZation 

HZSEL 2 r°; tek , es p,ace by U8ing a received si9nal of a known traini "9 sari - ■* 5! S3S 

40 a^uaT^^ 

^i^SSSL Q l er 21 fr ° m ° CCUmn9 in 3 proper manner 1716 accufa cy of the channel estimation can be 
2^ f ? f P ? P ^ rt,0n WhiCh th8 trainin9 Senes "^P* 3 °™ frame, but this degrades the trans 
befZnvSTr H ttle . , " tanded data - ^cording.* it is desirable that the accuracy of the channel estimation coufd 
be improved while reducint the proportion of the training series in one frame 

46 l° 010 L ^ fe " 0t ' imited t(> 3 reCeiver for multiple channel transmitted signals inclusive of MIMO but the same is 
thB^rf 16 ^ anne ' estima ''°n of a receiver such as RAKE receiver or a receiver using an adaptive aroy antenna where 
the certainty of a decoded result is improved by an iterative decoding process e array antenna wnere 

[001 1 ] The described turbo-receiver has following restrictions: 

50 ■ It is an accommodation for a single user (single transmitter) or only one series transmitted signal 

" acSairi. e t/„Sn a fln * iS , ne< f SSary in ^Producing an interference component, and this must be estimated in 
actual implementions. An estimation error results in a degradation in the effect of an iterative equalization. 

IS* J2 a " ° bjeCt I* 6 « V6nti0n t0 Pr ° Vide for these two restrictions by providing a turbo-reception 

55 of transmitttri • ""^ allow the receive ' ™"tioned a^ve to be extended to a receiver for a^raSy 

ron«, , T S ' 9na,S SUCh 38 f0r multiple usere or P aral,el transmissions from a single user 

He of a reSsiotatS 1?^?°" 5 ^ ! reCePti ° n meth0d and 3 r3ceiver "-^ in ««* a ^annel 
value of a receded s.gnal is estimated from the r ceiv d signal and a known signal s rving as a reference signal, the 
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receiv d signal is processed using the stimat d channel value, and the processed signal is decoded so that the 

al7ow r^ S,n9 . 7 Sti r e , d ? anne ' V8lU 8nd *° d6COdin9 are iterated U P°" *• -™ raceiv Tsjnal and whS 
allow the channel est.mat.on to b achi v d with g od accuracy using a relatively short known signal 

SUMMARY OF THE INVENTION 

KJL !2^fl?- to 8 -T 3SPeCt ° fth ° PreSent inVenti0n • there is provided a t"*o-reception method of receiving 
J£? h r" ,n , te9e T" al t0 or 9 reater ma " 2 > transmitted signals. The method comprises calculating a 
JrS ue >n(q) n= 1. N) from M received signals r m (m= 1, -. M) and N series known sfcnals. determinfng 

me d^o" ^" !k T b0 ' b * n(k) ° n thS ba8iS ° f N Series 8 P riori informati °" ^ n (k)l which are obtTed by 
me decoding, and us.ng the channel value h mn (q) and the soft decision transmitted symbol b' n k) to calculate an inte * 

ESSr^T !*' {k 2 f T ed by a " interSymbo1 interfa ^ Produced by an n-th tranl tfe^ £na iiel Tnd 
transmitted signals other than the n-th transmitted signal as follows: 



20 



H= 



H(o) ••• H (Q _ 0 
0 H (0) 



'(Q-l)J 



25 



30 



H (q) = 



l ll(q) 



h lN(q) 



_ h Ml(q) hMN(q) 



35 



B'(k)=[b* T (k+Q-1)...b' T (k)...b' T (k-Q+1)] T 



40 



b'(k+q)=[b , 1 (k+q)b' 2 (k+q)...b' N (k+q)] T 



q=Q-1 -Q+1 forq*0 



45 



50 



b'fkMb', (k)...O...b- N (k)] T q=o 

Tr b o\ h Z a ,7^T e T 8t r tH ^T' ° re P resente * enumb « r o^"'«Pathsofeach transmitted signal wave, 
m a L £ . m 1 re P re r Sents a tran *posed matrix, subtracting the intersymbol interference HB'(k) from the received 
matnxy(k) to obtain a difference matrix y'(k) where wreeeivea 

y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 



55 



rfrMr, (k)r 2 (k)...r M (k)] T 

enSSl 9 ^, 3 ? fi ' ter C ° effiCient Wn(K) aPP " ed l ° the received si 9 nal for the ""t" transmitted signal in order to 
tZT^m l '"France component in the difference matrix y'(k) using the channel matrix H or reference signaT 
filtenng th drtf rence matnx /(k) according to the adaptive filt r coefficient w n (k) to obtain a log-likelihood ratio for 
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bZEFUlT"* ds ' 9nalfr * en -* 8 ries transmitted signal from which th interference has been lim- 

inated. anddecoding by using the th log-lik lihood ratios for N series 

112 i S aSee ° ndaSPe,,,h Prasentinvntion.inth arrangement according to the first asp ct fo the 



b'(k)=[b' 1 (k)...-f(b l n (k))...b" N (k)] T 

« ?1* ki'T^^I" fe l0C3ted at P0SiU0n ' and * > resents a function which satisfies f(0) = 0 and having 
as a vanable b' n (k) which satisfies d{f(b'„(k))}/d{b , n (k)}s 0. 9 

[0016] According to a third aspect of the present invention, the equalization takes place in a plurality of staaes and 
£e number of senes in the equalization output is sequentially reduced in the later stages 

[0017] According to a fourth aspect of the present invention, in a turbo-reception method in which a channel value 

sLnaT^'f 'h IS ^ reC6iVed Si9nal and a *™ si9nal seeing as a reference sTgnaUne rocSeo 

s.gnal ,s processed us.ng the estimated channel value, performing a decoding processing upon the procesLSal 

^ ve^^ V8lUe ^ d " 0 *« ?£«—"» « «era te d X same' 

vaSHf an asso?l d 2ft h r hard deCiSi ° n infomia t ion si 9"al has is determined on the basis of the 

tetTua fto«^S2SS softdec,s,on ^fon signal, anda hard decision information signal having a certainty which 
.sequel to or greater than a given value Is also used asareference signal in the channel estimation o?the next iteration! 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0018] 

rSnLL 8 ^? ^T 9 , 1 tmataml arran ^nt <* * system including an embodiment of a turbo-receiver 
according to a first aspect of the present invention; waiver 

3? sh 2 own S?T? "° VfeW Sh ° Win9 8 SPe ° ifiC eXamP ' e 01 8 fUnCti0na ' arran9eman, of a mu, f P's output equalizer 
invention? ^ ° f ^ emb ° diment ° f a t"*o«ception method according to the first aspect of the present 
Fig. 4A is a diagram of an exemplary frame; 

^\^ZTl!!^T X l af ProceSsin9S wnicn ara formed during each iteration, presented in order to illustrate 
an .terabve channel est.mat.on method according to a fourth aspect of the present invention 

be^rlaiS; " eXemp ' ary fUnCtional """Cement for deriving a hard decision symbol which is likely to 

Fig. 6 is a flow chart of an exemplary processing procedure for channel estimation according to the present inven- 

Fig. 7A is a diagram of an exemplary functional arrangement of part of the equalizer 31 according to the second 

aspect of the .nvention which reflects an error correcting decoding result of a signal being £2d- 

Fig. 7B is a diagram showing an exemplary processing procedure therefor; 

Fig. 8 is a diagram of an example of a receiver which uses a turbo-equalizer iteratively 

Fio foV.'lEZ.l- exemplary receiver which performs RAKE reception-turbo-decoding iteratively; 

feraJve" 6XemP * rece,Ver wnich P erforms a " «*P*" array antenna reception-turtJ-decoding 

Fig. 11 A is a schematic view of a turbo-equalizer; 

Fig. 11 B is a schematic view of a turbo-decoder; 

Fig. 12 is a schematic view of a receiver which iterates a processing of a received signal using an estimated 
channel and a decoding processing of the processed signal; esomatea 

?« I 3 ^H° W Ch , art ° f a " exemp,ary Passing procedure of a reception method which iterates a processing 
of a received s gnal usmg an estimated channel and a decoding processing of the processed signal; 9 
Fig. 14A is a diagram of an exemplary frame arrangement; «»wnw. 

S«» f. SL^IT"? i, ' US ! rati ° n of an iterative Processing of an estimation of a channel H and a noise covariance 
matnx U when a received signal contains noises other than white Gaussian noise- 

es«m 1 a 5 «on a o f SfT 34 ' 0 Vi8W ° f a " exemp,ar y functional arrangement of part of an'equalizer which is used in the 
estimation of the noise covanance matrix U; 

inlhlestiml 0 ^^ a " eXemp,a,y Passing procedure which iterates the estimation of a channel value used 
in the est.mat.on of the noise covanance matrix U and the decoding processing- 

Fig. 17,s a diagram showing th principl of a turbo-r ceiv r according to a third asp ct of th pres ntinv ntion- 
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Se preset REif ^ XemP ' ary fUn ° ti0nal arran9ement of 8 turbo-receiver according to th third aspect of 
^s^;^ r(prestage) qualiz r 

2£5 Z^&£r*" Pr0CeSSin9 PrOCedUre ° f 3 ^»P«- -thod acceding to the third 

SoSi £2K ipr^^ afran9ement of a mM9e e9Ualizafion - cUon 

S 9 thf„~L di f ? ram r f a ?o e . XamP '^ ° f 8 SySt6m arran 9 eme »» which an embodiment according to the first aspect 
of the present invention (2) is applied; ^ 

SSL22^ ka S r Sh ° WS «" 6 T rate charactori8,lc °f a turbo-receiver according to the first aspect of the present 
7Z2Z ^ °" ^ Ch8nnel ^ b6en PerfeCUy eStimated ' ^ reDresenti "9 a "it powerand N 0 

characterisec when - itera,ive *— esumauo - has >*«™<* 

Fig. 25 graphically shows an error rate characteristic of a turbo-receiver according to the fourth aspect of the 
present invention, in particular, employing the iterative channel estimation- 

colariancTmaMx U; ° WS " ^ "** CharaCteriStic of the ^receiver which employs an estimation of a noise 
Fig. 27 graphically shows an error rate characteristic of the turbo-receiver shown in Fig 1 • 
S"^L B H n 2I C !!! y f?K WS a " eTOr ^ characteristic as P'°tted against N 0 of an embodiment according to 
PteSe °" WNCh refleCtS a " erTOr COrrecti " 9 decodin9 result of a si 9" a ' being 

2Z^^Zl^ullT Ulati ° n ° f a " en0r ^ cha — -fa turto-receiver according to 
Fig. 30A is a diagram illustrating the concept of a MIMO system- 
ic? 1 r dia9ram °? n arran9ement in which received signals from a pair of antennas are input to a turbo- 
receiver as four received signal series; and 

Fig. 31 is a diagram of a functional arrangement of a conventional turbo-transmitter and receiver for a single user. 
EMBODIMENTS OF THE INVENTION 
35 First aspect of the Invention ( 1) 

2 ^JcS^ exemplary arrangement of an MIMO system to which the present invention is applied. 
2S2.J . transm Jtters 31 -SN. information series c,(i) to c N (i) are encoded in encoders 11-1. ••• 11-N 

ttus modulating a earner signal in accordance with these modulation signals to transmit signals b l( k) to b N (k) In thfs 
SET STl E 8,9na 1 ^: (k) : bN(k) fr ° m the «. ™.SN form N series transmitted signal * 

L fnout to a ^ f "T ^ W , hfeh iS reC6iVed by 3 mU ' tiple 0Utput receiver mrou 9 h trensmission paths (channels) 
s Zn i°l?? P t °Tl I qUa,,Zer 31 * A Si9na ' reCeived by 018 receiver is «««« in »o a baseband signal, which 

^o Q r m p , efe afe 006 ° r more di9ital si9nals - 016 numt >er of which is represented by an 

signals " " '* ^ M *"> f0fTned in, ° received s, ' 9nals in *. form of M digital 

[0022] The equalizer 31 delivers N log-likelihood ratios A, [b,(k)J. ■•• A, [b N (k)l In subtracters 22-1 - 22 n a nri^n 

whllh:! w ' nf ^ l f SOf * ,npUt SOft output < S,SO > dec oders (channel decoders) 24-1. ■». 24-N to be decoded 

vers 26-1 -, 26-N to 
, 22-N, respectively. 
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rmM-Icpo 0-1 2n=i N h mn (q)b n (k-q) + v m (k) 



(20) 



where q-0,- Q-1 , and Q represents the number of multipaths for each transmitted wave. Defining a matrix y(k) by a 
5 similar procedure as used for the single user, we have 



10 



15 



where 



y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 
=HB(k)+n(k) 

r(k)=[r 1 (k)...r M (k)] T 



(21) 
(22) 



20 



25 



where 



H= 



H(o) 
o 



(Q-l) 

H(0) 



'(Q-l). 



(23) 



30 



H(q) = 



n H(q) 
h Ml (q) 



h lN(q) 
h MN(q) 



(24) 



40 



45 



SO 



55 



B(k)=[b T (k+Q-1 )...b T (k)...b T (k-Q+1 )J T 

b(k+q)=[b 1 (k+q)b 2 (k+q)...b N (k+q)] T 
q=Q-1,Q-2 -Q+1 



(25) 



(26) 



mSS m th ln,erfe ; ence el ' m '"ation step, it is assumed that it is desirable to obtain a signal from an n-th user (trans- 
m.tter). n his example, a soft deas,on symbol estimate for signals from all users (transmitters) and a channel matrix 

EnfZTIn'* ,mP :' S ? T* 0 ™ matriX) " are US6d 10 produce a *^°* s of interferences by user signaToZ 
2£! ^ . ^ 8er f- lnterferences created "y «ser signal itself, or an interference replica H-B'(k) is repro- 

duced. Then, the interference replica is subtracted from y(k) to produce a difference matrix y'(k). 



y'(k)=y(k)-HB'(k) 



=H(B(k)-B'(k))+n(k) 



(27) 



(28) 



where 
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B , (k)=[b ,T (k+Q-1)...b ,T (k)...b' T (k-Q+1)] T (29) 

b , (k+q)=[b , 1 (k+q)b , 2 (k+q)...b , N (k+q)] T 

:t l = Q-1 -Q+Lq^O (30) 

b , (k)=[b' 1 (k)...O...b , N (k)] T : q=0 {31) 

SmS wh^ tT"' 3 I k * T'f ?■ K '! t0 be understood *at b' n (k) represents a soft decision transmitted symbol 
me Inherent, symS " ^ ^ maWx B ' (k) represents a rep,ica ma « x of 

namfiv ^.JE.!?!^ ""^l™ *? WhiCh * USed to eliminate me residue of *° interference component, 

namely, a res.dual interference based on the imperfectness of the interference component replica H-B'(k) and an in- 
terference component created by the n-th signal itself .is determined according to MMSE (minimum mean square enTr) 
catena as one which minimizes the following equation (32): ' 

w n (k)=arg min||w n H (k)y(k)-b n (k)|| 2 {32) 

[0026] Subsequent operation remains the same as for the single user. Specifically. w n (k) which is obtained in this 

l"T 1°* I"' 318 r nH(k) y " (k) ' and a resu,t ° f is fed through deleaver 23-n to be Spu as X? 

[b„(i)] to the decoder 24-n where a decoding calculation is made. 1 
[0027] The described method of applying a filter (linear equalization) processing upon the received signal r is re- 
aH tirasA^r) rtn ^ u ®ers ^ ^h^h 8 consequence, the number of outputs from the equalizer 31 will be e^uarlo'N, and 
a I these outpu s are decoded by respective decoders 24-1 . 24-N. What has been mentioned above is an extension 
of a turbo-receiver for single user to a receiver for multiple users (MIMO) 

%Szj£Zm£l7T 9 deS c crip " 0 "' " w "' beseen * aXan exemplary functional arrangement of the multiple output 

matn^H whi^lT, h ,S | ? enerated and is fed to e 1 ualiz «™ 312-1 to 312-N for each user. The channel 

EELS I J ,^ ^ 3 ° hannel eSt,mat0r 28 is a,so su PP ,fed to the equalizers 312-1 to 312-N. A priori 

222 ?- ^ W K 6 f Ch ° hannel deC ° der 24 " n i8 SUppHed t0 a ~ n decision ^bol estimator 313 where a soft 
decision ttansm.tted symbol esfmate b' n (k)= tan h [b n (k)]72] is calculated. All of equalizers 312-1 to 312-N have 

will b^Sed arran9ementS and process in an identical manner . end accordingly, a typical one of them (312-1) 

mlTJn^Sr ° 1 n l S i° ft dedSi0n ^ smitted s y mb ° l esfimat es *,W to b' N (k) are supplied to an interference replica 
2 m T R Lv r 3 " lnterference re P |ica ma tnx B- l( k) is generated according to the equations (2 9 ) to (31). 
tan the matnx B l( k) ,s subject to a filter processing according to the channel matrix H in a filter processor 315-1. and 

316-1 to produce a difference matrix y'^k). »»w«ih 
SSL w^"! n C ! ,anne ' maWX H ° r 8 referenc e signal which will be described later is input to a filter coefficient 

ofnTnt 1 tl to de ^ fi ' ter C0effiCi6nt Wl(k) Which iS USed to aliminate 016 resid "e of the interference com- 

ponent. In the example shown, the channel matrix H and a covariance a* of a noise component and soft dedsion 
transmitted symbols b' l( k) to b' N (k) from the soft decision symbol generator 313-1 ^re^^J^^S 

h£™J? Z . ?" PeCrfiC eXample ° f d etermining the filter coefficient w,(k) will be described later. The 

difference matnx y l( k) ,s processed with the filter coefficient w l( k) in an adaptive filter processor 31 8-1 and A [b Z 
,sde nrered as an equa.ized output for the received signal which corresponds to the transmitted signal frem userl . 
J I L pr ° cess ' n9 P roce d"re for the multiple input multiple output turbo-reception method according to the de- 
scribed embodiment of the present invention is shown in Fig. 3. At step S1. a channel value h mn (q) and a covariance 
chateiZ^T" a ™, Calcu ! ated a receiv «f signal r(k) and each training signarG). At step S2. the 

PSTJSS? ! *££ T, ?H Cha r ' V K a ' Ue hmn(q) - A< SteP 831 016 80,1 decision transm '« ed symbol estimate 
b „(k)- tan h ^ [b„(k)j72] is calculated on the basis of the a priori information ^ [b„(k)] which is obtained during a 
previous iteration of the turbo-reception processing. " 9 

[0032] At step S4. the received signal y (k) is general d from the receiv d signal r(k). At st p S5. th int rference 
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replica matrix B' n (k) is general d according to the quations(29)to(31)usingth respective soft decision transmitted 
symbol stimat s b„ (k) . At st p S6. the int rf rence component replica H-B' n (k) for th receiv d signal frem th n- 
« nS T^ ' S calculated - At st P S7 > *e interference component replica H-B' n (k) is subtracted from th receiv d 
matrix y (k) to provide , the diff rence matrix y n (k). At step S8. using the channel matrix H. the soft decision transmitted 
symbols bj(k) to b N (k) and the covariance o* of the noise component, the filter coefficient w n (k) which is used to 
eliminate the residue of interference remaining in the received signal from the n-th transmitter is determined according 
to the minimum mean square error criteria so as to minimize the equation (32) 

5°!? ■ » T P ?l ■r? 6r processin 9 acco "*n9 to the filter coefficient w n (k) is applied to the difference matrix y n (k) 
to obtain the og-l.kel.hood ratio A 1 [b n (k)]. At step S10.. the a priori information A, [b n (k)J is subtracted from A, [b"(k) 
and the result is demteneaved and decoded to deliver the log-likelihood ratio A 2 [b n (k)J. The steps S4 to S10 are 
Performed either simultaneously or sequentially for n=1 to N. Subsequently, at step S11, an examination is made to 
see whether or not the number of decoding operations, namely, the number of turbo-reception processings has reached 
a given numberlf the given number is not reached, the operation goes to step S12 where the extrinsic information X, 
[b n (k)] is subbacted from the log-likelihood ratio A 2 [b n (k)], and its result is interleaved to determine the a priori infor- 
mation ^ [b n (k)]. thus returning to step S3. If it is found at step S1 1 that the number of decoding operations has reached 
the given number, a resulting decoding output is delivered at step S13. 

[0034] The channel estimator 28 will now be described. Each received signal r m (k) can be represented as follows: 

r m (k)=5: q= o Q ' 1 £n=i N h mn (q)-b n (k-q)+v m (k) (33) 

[0035] The channel estimator 28 determines the channel value (transmission path impulse response) h mn (q) appear- 

i?kll?r <33) *! T" P ° Wer (Sa2) ° f the n ° iSe V -" (k) " Norma "* a uniaua word (fining s?gnal) P which 
is known to the receiver is inserted at the beginning of each frame to be transmitted on the transmitting side, as shown 
in Fig 4A. and the receiver utilizes the unique word (known signal) as a training series to estimate the channel value 
n mn (q) using RLS (recursive least square) technique. Each of the channel decoder 24-1 , -, 24-N delivers +1 if a log- 
likel.hood ratio A 2 [b, (i)]. • A 2 [b N (i)] ispositive* and -1 if the latter is negative as a decoded code signal (or transmitted 
encoded symbol hard decision value) b,©.™, b N (i). which are fed iteratively through interieavere 27-1 ••• 27-N to the 
channel estimator 28. The received signal r(k) is input to the channel estimator 28 together with a unique word which 
is fed from a unique word memory 29 to serve as a reference signal. On the basis of these input signals, the channel 
estimator 28 estimate each h mn (q) according to the equation (33) and o* value according to the least square technique 
This estimation can be made in a similar manner as the estimation of an impulse response when adaptively equalizing 
JnA,^ S ' 9na ' Wth a " ada P ttve filter bv estimating an impulse response of a transmission path 
[0036] The use of the training series in this manner is conventional, but in oraer to enhance the net transmission 
rate, it is necessary to reduce the proportion that the unique word occupies in one frame, but this increases the error 
of the channel estimation. If there is such an error, it results in degrading the iterative equalization response mentioned 
posed 9 prBSent inventi °n. an iterative estimation of the channel value in the following manner is pro- 

[0037] The concept according to the present invention is illustrated in Fig. 4B. The purpose of the conceDt is to 
estimate toe channel value iteratively in each stage of the iterated equalization processing of the same received signal 
or dunng the each iteration of the turbo-reception processing. During the first run. for information symbol series which 
follows the unique word, the channel value is estimated using only the unique word as a reference signal, and the 
estimated channel value is used to equalize the received signal and to estimate the transmitted symbol However 
before the second iteration of the equalization processing, the channel estimation is made using the unique word as 
the reference signal, and the symbol estimate (hard decision value) which is obtained during the previous decoding 
processing is also used as the reference signal to perform a channel estimation of the entire frame. In this instance 
rather than using every hard decision value, only that hard decision value or values which have been determined as 
■ke y to be certain is used as the reference signal. The hard decision is rendered by delivering +1 when the logic- 
l.kel.hood ratio A 2 [b n (.)J from the decoder 24-n is positive, and delivering -1 if the later is negative. The greater toe 
absolute magnitude of the log-likelihood ratio A 2 [b n (i)], it can be said the greater the likelihood that the hard decision 
value is likely to be certain. For example, the certainty that is rendered when the log-likelihood 5 is determined to be 
+1 is greater than the certainty when the logic-likelihood 0.3 is determined to be +1 . Accordingly, an interative channel 
be deS ribed tetow US6S 8 threSh0W Va ' Ue to desi 9 nate h ard decision value b n (.) which is likely to be certain will 

[0038] Initially, using the log-likelihood ration A 2 [b n (i)] from the decoder 24-n. a soft decision symbol value b' (I) is 
determined as follows: nv/ 
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b' n (j)=tan h[A 2 [b n <i)]/2] 

This operation is made in ord r to normaliz th logic-likelihood value to 1, thus preventing its absolute magnitude 
from exceeding 1 . Next, a threshold between 0 and 1 is provided, and when the absolute magnitude of a soft decision 
value b n (i) is greater than the threshold value, a corresponding hard decision value b n (i) is preserved to be used in 
the iterative channel estimation. For example, if the threshold value is chosen to be 0.9. those of hard decision values 
b n (i) which correspond to soft decision values b' n (0 having absolute magnitudes equal to or greater than 0.9 are se- 
lected. It is believed that the certainty of the selected hard decision value b^i) is high in view of the fact that the threshold 
value is as high as 0.9, and accordingly, it is believed that the accuracy of the iterative channel estimation which is 
made using such hard decision values can be improved. However, it is also considered that because the number of 
selected symbols is reduced, the accuracy of the iterative channel estimation may be influenced and becomes degrad- 
ed. Accordingly, it is necessary that an optimum threshold vajue be choosed between 0 and 1. To add. if the threshold 
value is to be 1 . if follows that there is no hard decision value b n (i) which is selected, stating to the effect that no iterative 
channel estimation takes place. Accordingly, a threshold value on the order of 0.2 to 0.8 is chosen, as will be further 
described later. 

[0039] Accordingly, those of transmitted symbol estimates (hard decision values) Mi),»,b N (i) for the information 
symbol senes during the first transmission which are determined to be likely to be certain according to the threshold 
value are fed from the outputs of the interleaves 27-1, 27-N to a previous symbol memory 32 to be stored therein 
as a previous transmitted symbol estimate. During the second iteration of the equalization and decoding processing 
for the received signal r(k) (it being noted that the received signal r(k) is stored in a memory), the unique wo,rd is initially 
used to make a channel estimation, and those of estimated transmitted symbol hard decision estimates b«(i) ■■ b N (i) 
which have been determined to be likely to be certain are read from the previous symbol memory 32 to be fed to the 
channel estimator 28 to make a channel estimation, namely, a channel estimation for the entire frame A resulting 
estimate h mn (q) and a 2 are used to perform an equalization and decoding (the estimation of transmitted symbol) with 
respect to the received signal r(k). At this time, those symbol values among the estimated transmitted symbols which 
had been determined to be likely to certain according to the threshold are iused to update the content of the previous 
symbol memory 32. Subsequently, an estimation using the unique word and an estimation using those of the previously 
estimated transmitted signals which are determined to be likely to be certain are used to make a channel estimation 
of the entire frame during a subsequent iteration of the equalization and decoding. The estimated channel is used to 
perform the equalization and decoding (the estimation of transmitted symboB and to^pdate the previous symbol mem- 
ory 32. Alternatively, those of the transmitted symbol hard decision values b,(i).- t b N (i) from the decoders which are 
determined to be likely to be certain according to the threshold value may be directly stored in the previous symbol 
memory 32 to update it, and when the symbol values stored in the previous symbol memory 32 are to be utilized they 
will be fed through the interleaves 27-1 . 27-8 to be input to the channel estimator 28. 

[0040] By the iterations which proceed in this manner, an error of the channel estimation can be reduced the accuracy 
of the symbol estimation can be improved and the problem of a degraded response in the turbo-equalization due to 
an error of an channel estimation can be improved. 

[0041] When a channel estimation is made with respect to information symbol series using the symbol hard decision 
values which are likely to be certain in the manner mentioned above, a function as indicated in Fig. 5 is added to each 
decoder 24-n. The log-likelihood ratio A 2 [b n (i)] is input to a soft decision value estimator 241 to calculate b'(i)= tan h 
[A 2 [b n (i)]/2] thus estimating a transmitted symbol soft decision value b' n (i). The value b' n (i) is compared against a 
threshold value Th from a threshold presetter 243, thus delivering 1 when b' n (i) is equal to or greater than Th and 
dehvenng 0 when it is less than Th. On the other hand, the log-likelihood ratio A 2 [b n (i)] is input to a hard decision unit 
344 which delivers +1 if A 2 [b n (i)] is positive and delivers -1 if the latter is negative as a symbol hard decision value 
b n (i). The symbol hard decision value b^i) is fed through a gate 245 which is opened when a corresponding symbol 
soft decision value is equal to or greater than the threshold value, and is passed through the interieaver 27-n shown 
in Fig. 1 to be fed to the previous symbol memory 32, thus updating the stored symbols. 

[0042] A procedure of the channel estimation which also utilizes a symbol hard decision value or values which are 
likely to be certain is shown in Fig. 6. Initially, at step S1 . a channel estimation is made with a received signal r(k) and 
a unique word. At step S2. an examination is made to see if the decoding processing is for the first time, and if it is at 
step S3, the estimated channel value h mn (q) is used to perform the equalization and decoding processing or the od- 
erations shown at steps S3 to S1 0 in Fig. 3. 

[0043] At step S4, a transmittecteymbol hard decision is rendered with respect to the log-likelihood ratio A, [b (i)l to 

? t"?-!? 6 ^ h3rd dedSi0n V3,Ue bn(i) * At Step S5 ' b 'n (i)= ten h ' A 2 fenW 2 ! is calculated from the log-likelihood ratio 
A 2 Jb (.)], thus estimating a transmitted symbol soft decision value b' n (i). At step S6, by examining whether the symbol 
soft decision value b' n (i) is equal to or greater than the threshold value Th or not, those of the symbol hard d cision 
valu s b n (i) which are likely to be certain are d t rmined. At step S7, the symbol hard decision valu or values which 
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are likely to be certain are us d to update th cont nt of th previous symbol memory 32. At step S8, an examination 
is made to se whether or not the numb r of decoding op rations has reached a giv nvalu , and if not. the operation 
returns to step S1 , or more exactly, going through th step S1 2 shown in Fig. 3 and returning to the step S1 shown in 
Fig. 3. 

[0044] If it is found at step S2 that the decoding processing is not for the first time, at step S9, a previous stored 
symbol, namely, a hard decision symbol which is likely to be certain is read from the previous symbol memory 32 and 
it is used together with information symbol series of the received signal r(k) to make a channel estimation thus frans- 
femng to step S3. 

[0045] In the above description, even during a second and a subsequent iteration, the channel estimation is made 
from the beginning using the unique word as a reference signal, but it is also possible that only the hard decision symbol 
which is likely to be certain may be used as the reference signal during a second and a subsequent iteration In this 
instance, as indicated in broken lines in Fig. 6, an examination is made to see if this processing is for the first time at 
step SV, and if it is, the unique word is used as a reference signal together with the received signal to estimate the 
channel value. After storing the estimated channel value and parameter values which are used in the estimation in a 
HEmS? 31 StGP S3 '' the operation tonsfers to *e equalization and the decoding processing which take place at step S3 
[0046] If it is found at step SV that the processing is not for the first time, before the channel estimation is made a 
previously stored channel estimate and various processing parameters are set up at S4\ whereupon the operation 
transfers to step S9. 

[0047] It is to be noted that a solution of the equation (32) is given as follows: 

w n (k)=(HG(k)H H +o 2 I)" 1 h (34) 

where I represents a unit matrix, o 2 an internal noise power of a receiver (covariance of noise component) o 2 I a 
covanance matnx of noise component and G(k) corresponds a squared error of the channel estimation. 

G(k) =E[(B(k) -B'(k) )-(B(k) -B'(k) ) H J 

=diag[D(k+Q.l),. . . J}(k) 9 . . .,D(k-Q+l)] (35) 

where E [ ] represents a mean, diag a diagonal matrix (having elements which are not along the diagonal being zero's). 

D(k+q)=diag[1 -b' 2 1 (k+q) 1 -b' 2 n (k+b) 1 -b* u (k+q)] (36 ) 

q=Q-1, Q-2 -Q+1, q*0 

and when q=0, 

D(k)=diag[1-b' 2 1 (k) 1 1-b* N (k)] (37) 

T appearing in the matrix D(k) represents an n-th element (it is assumed that an n-th user's transmitted siqnal is 
assumed to be a desired signal). 
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h = 



H 



1 . (Q-l) • N+n 



H 



2.(Q-l)N+n 



L^M.Q.(Q-l).N+n 



(38) 



[0048] Thus, h comprises all the elements in the (Q-1)N+n-th column of H appearing in the equation (23). Input to 
the nuse power o* which are estimated in the channel estimator 28 and the soft decision transmitted symbols b\(k) 
SKcS^ 

w n (k) according to the equations (34) to (38). 

[0049] It will be seen that the equation (34) requires an inverse matrix operation, but the required calculation can be 
reduced by using the matrix inversion lemma for the inverse matrix. Specifically, by approximating each b* appearing 
in the equations (36) and (37) by 1. it follows that 9 appeanng 



D(k+q)=diag[0 0]=0 (q*0) 



(39) 



D(k)=diag[0 1 0] 



(40) 



SSL. 1"*' ? W h ! Ve e ' ement , S haVi " 9 3 ValUS ° f 1 for those elements are located at ™» ^d n-column 
while all ofter elemente are equal to 0. When the error matrix 6(k) of the equation (35) which is determined by the 
equations (39) and (40) is substituted into the equation (34), we have 



w n (k)=(hh H +o 2 l)- 1 h 



(41) 



where h is as defined by the equation (38). 

[0051] With this approximation, w n (k) does not depend on k. and accordingly, an inverse matrix calculation at every 
discrete bme k can be dispensed with, thus reducing the quantity of calculation 

[0052] The matrix inversion lemma for the inverse matrix is applied to the equation (41). The lemma states that 

^THX»ZZ M) T Tf ^ ° 8 (M> N) maWX and ° 3 (N " N) ^ maWx ' A SEE 
A '-B + CD-^", the inverse matnx of A is given as follows: 



A' 1 =B-BC(D+C H BBC)' 1 C H B 
Applying the lemma to the inverse matrix operation appearing in the equation (41). we have 

h(k) h(k) H +a 2 l=B' 1 +CD" 1 C H 



(42) 



h(k)h(k) h 



=CD- 1 C H .o 2 l=B 1 ,h(k)=C 



=D" 1 ,h(k) H =C H 



Using these equations to calculate the equation (42). the inverse matrix operation appearing in the equation (41) can 
be calculated. While the equation (42) contains an inverse matrix operation (D + CH B BC)-1 . this inverse matrix becomS 
scalar and can be simply calculated. 
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[0053] Thus, in this instance, it is reduced to th following form: 

w n (k)=1/(o 2 +h H h)h (4^!) 

1/( ) on the right side of this equation is scalar or assumes a constant value, which may be chosen to be 1 Thus we 
can put w n (k) -h. indicating that w (k) can be determined by only h. As indicated in broken lines in Fig. 2, only h which 

c^Sm^matciau-T 00 " (38) the Channe ' maMx H * om * 9 channel estimator 28 ma * be * ■» 

rtlLS^ 6 '?, 9 ^^ 10 USe ° f the matdX inVerSi ° n lemma for the inverse maWx for me approximation by 
fte equabons (39) and (40). the approx.mation also allows the quantity of calculation for the equation (34) to be reduced 
In pamcular when th.s approximation is used, and the matrix inversion lemma for the inverse matrix is used, the quantity 
ZJSmZT ^ « fUrth8 :rl! JC8d ' additi ° n ' ifthB °° variance matrix ° f noise component is represented 

l iZ^ « 1 W S (k) =h be " Sed 38 indi0ated by 1,16 equation f 41 " 1 )- whereupon it is independent 

from the covanance matnx. allowing a further simplification of the calculation. 

Second a spect of the invention (reflecting an error correction) 

[0055] In the equalization processing where H-B'(k) is subtracted from the received matrix y(k) represented by the 

SlSlfS a " e r^ eCti °? deC ° ding r6SUlt iS reflected in a ^nsmitted signal soft decision value for a signal 
nlwt h f Snal b„(k) be,ng detected, but an error correction decoding result which relates to the signal b n (k) being 
detected is not. For this reason, it is preferred to employ the following processing ' * 

[0056] b'(k) appearing in the equation (29) or the equation (31) is changed as follows; 

b'M-PtWb'iM... 
b'„.i (kH(b' n (k)) b' n+ i (k)...b' N (k)] (43) 
where f(b* n (k)) is an arbitrary function having b' n (k) as an input. 

S^JSJT^ • ChB ^Ti * beC ° mes P 088 "" 6 10 reflect an error correction decodi "9 result with respect 
to the signal b n (k) being detected. Thus, rather than using b' n (k)=0. by adding a suitable value depending on b'(k) it 

to KX'SSS* 3 Si9na ' bSin9 deteCted Wh ' Ch iS bUned in 3 n ° iSe ° f a " interferin9 8i9nal ' thUS a ' lowi "9 b n (k) 

t«hwL B T USe S l 9 K n ? b * n(k) relateS 10 3 reSU " 0f a nard decision rendered "P° n a symbol which corresponds 
£££ S^T k ^ 816 9feater the abS ° IUte ma9nitude of b '" (k >' me hi 9 her tha reliaW «ty of the hard 

JUT™, Tr * , ~ rTesp< : nd8 40 b ' n(k) ' * iS neCeSSary that f(b « (k » saUsfies 1,16 followi "9 requirements: 
[0059] The function f has a value of 0 for b' n (k)=0 or when the reliability of the hard decision symbol is equal to 0. Or. 

f(0)=0 (44) 
[0060] In addition, the greater the value of b- n (k). the greater the value of the function f. Or, 

<Kf(b\,(k)))/d{b' n (k)}*0 (45) 
Examples of such f(b' n (k)) include the following: 

f(b' n (k))=<xxb' n (k) (46) 

f(b , n (k))=axb' n (k) 2 (47) 

^rit!^ eqUati °" ( . 46) / S . USed and a iS Ch ° Sen t0 a COnstant - the e< > uation < 43 > «" be implemented 
in a simple manner. Here, a must satisfy the requirem nt: 0<o<0.6. If a is greater than 0.6. BER (error rat ) will b 
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d graded, prev nting a correct decoding result from being obtained .It is also contemplated to mak a tab variable 
.naccordanc with the reliability of a d coding result. For xampl . a may be chosen for achit ration of the decoding 
^T 9 '!" tn,s,nstence - OTe Pliability of the d coding result will be normally improv dwiththenumb r of iterations 

^2« n n! C °rr 9 H ProC ! SS,n9S f° COrdin9ly - an increased value mav be for a depending on the number of 

iterations of the decoding processing. Alternatively, the reliability of an entire decoded frame may be determined upon 
each iteration of the decoding processing, and the value of a may be chosen in accordance with the reliability thus 
SlTwh h I emi " e , the re,iabilitv of the decoded «*me, a decoding result may be compared against a decoding 
result which is obtained dunng an immediately preceding iteration, and a number of hard decision symbols which 
changed from the previous decoding operation may be counted, for example. Thus, where there is an increased number 
of changed decision symbols, the reliability may be determined to low while when the number of changed hard decision 
symbols is small, the reliability may be determined to be high. 

Ew 1I *««=? (l ? ' S Cha " 9ed in thte manner ' the equation (35) whicn is used when determining the filter coefficient w„ 
(k) for MMSE (minimum mean square error) filter is desirably changed as follows: 

G(k)=E[(B(k)-B'(k))(B(k>B'(k)) H ] 

=diag[D(k+Q-l),..., D(k),..., D(k-Q+1)] 

[0062] Using the equations (2g) and (31 ), it follows that assuming 

B'(k)=[b'(k+Q-1)b , (k+Q-2) ... b'(k) ... b'(k-Q+1)] T 

b'(k)=[b' 1 (k)b l 2 (k) ... -f(b' n (k)) ... b' N (k)] T 

b , (k+q)=[b' 1 (k+q)b' 2 (k+q) ... -f(b' n (k+q))... b' N (k+q)] T 
:q^0, q=Q-1 -Q+1 

D(k) has elements located at n-row and n-column, which are represented as follows: 

E[(b„(k) + f(b' n (k)))-(b n (k) + f(b' n (k)))*J 
where ( ) represents a complex conjugate. For BPSK modulation, this expression turns into the following expression: 

E[b n (k) 2 +2b n (k)f(b' n (k))+f(b' n (k)) 2 ] 
=E[b n 2 (k)] + 2E(b n (k)f(b' n (k))+E[f(b' n (k) 2 ) 
The first term has a mean value of 1 . When b n (k) is approximated by b'(k). the equation (37) turns into the following form: 

Dfk^diagM-b-Vk) 1-b l2 2 (k)...1-b' 2 lv1 (k) 1+2E[f(b'„(k)b' n (k)) 

+ E[f(b' n (k) 2 ) 1-b ,2 n+1 (k)...1-b' 2 1 (k)] (4S) 
[0063] For example, when the equation (46) is chosen for f(b' n (k)). D(k) turns into the following form: 

D(k)=diag[1-b' 2 1 (k) 1-b ,2 2 (k)...1-b l2 n . 1 (k) 
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1+<2a+aV„W 1-b- 2 n+1 (k)...1-b ,2 1 (k)] {49) 

[0064] An exemplary functional arrangement which estimates an adaptive filter coefficient w n (k) when reflecting an 
error correct.on decod.ng result into a signal being detected is shown in Fig. 7A where a signal being detected is chosen 
as a fransm.tted s.gnal b,(k) from a first transmitter. A soft decision transmitted symbol b'^k) is input to a function 
calcu ator 331-1 to calculate a function f(b- l( k)). Soft decision transmitted signal b' l( k) to b' n (k) from N decoders ano f 
(b ,(k)) are input to an error matrix generator 332-1 to calculate and generate an error matrix G (k) according to the 
equabons (35). (36) and (48). The error matrix G(k), an estimated channel matrix H and noise power o* are input to a 
filtercoeffiaentgenerator333-1 where a calculation of theequation (34) is made to estimate an adaptive filter coefficient 
w (k). In this .nstance. f(b n (k)) is also input to an interference replica matrix generator 314-1. thus generating an 
interference replica matrix B'(k) represented by the equation (29) according to the equations (30) and (43). The filter 

^A° le nl JSPt"^ i 1 f !l t L rinS ? 3 difference matrix VW in ^ adaptive filter 318-1. thus yielding a log-likelihood 
rabo A, [b,(k)] It w.ll be noted that ,n the filter coefficient estimator 317-1 shown in Fig. 2. the function calculator 331-1 
shown in F'9- 7A 's omitted, and only the soft decision transmitted symbols b-,(k) to b' N (k) are input to the error matrix 
generator 332-1 in order to calculate the equation (34). 

OSL sVt? JTiS?/ R9 m tt,e J n,e ^ renCe reP ' iCa maWX B ' (k) iS 9enerated at step and after ^cessings 
at steps S5 to S7. the filter coeffiaent w n (k) is determined at step S8. If a calculation of the equation (34) is made 

EUSST^SI! 9 31 S ? P S8, S ° ft de ° iSi0n transmitted svmb °'s »>i(k) to b' N (k) are used to calculate the equa- 
KlSK ( ^ 9ene f a " erTOr maWX G(k) 31 Step S8 - 2 ' and * e erTOr matrix G < k >' estimation channel matrix 

In L ^ I P0Wef t a Z USed t0 Ca,CUlate the eqUation (34) to d6t «""ine an adaptive filter coefficient w n (k) at 
siep 08-3, as shown in Fig. 7B. n 

[0066] When it is desired to reflect an error correction decoding result in a signal being detected as mentioned pre- 
viously, a soft decision transmitted symbol b' n (k)of a signal which is to be detected may be calculated at step S8-1 
before entenng the step S4. and this may be used at step S4 where the equation (31) is replaced by the equation (43) 
or in other words, the equations (29). (30) and (43) may be used to generate an interference replica matrix B' (k) and 

S T^f (37) ma/ 06 rep ' aCed by the e(,Uation (48 >- ln the event *'nW) is chosen to be equal to a 

b n (k) or a b n (k)2 and when a is chosen to be variable, a may be determined in accordance with the number of process- 

%£^JnZg!KT Sntife deCOded *— 81 SteP S8 " 1 - 1 - and 1+ ^)b' n (k, 2 may be calculated 

EHf? * Th6 techniq " e °J u reflecti "9 an error correcting result into a signal being detected is also applicable to a single 
user turbo-rece.ver wh,ch has been described initially in connection with the prior art. In the technique of reflecting an 
error correcting result into a signal being detected, the approximation represented by the equations (39) and (40) can 
applied, and hi .this instance, only a matrix h shown by the equation (38) which is supplied from the channel estimator 
28 may be input to the filter coefficient generator 333-1 . as indicated in broken lines in Fig 7A 
[0068] In the above description, the adaptive filter coefficient w n (k) is determined according to the equation (34) or 
by using the channel matrix H. but the use for the channel matrix H can be dispensed with. Specifically, during the 
llnfllw I" 9 P ^ CeS . S,n i (t "*°- rece P tion Processing), the error matrix G appearing in the equation (34) becomes a 
dSln^nc » h 9 V ' ? ! d ? e : enCe . matrix *'< k > and the traini "9 ^9nal either alone or in combination with a hard 
decs on transmitted symbol b„(k). preferably b n (k) having a high reliability in the sense mentioned above are input to 
Uie -filter •coeffic.en generator 333-1to calculate the adaptive filter coefficient w n (k) in a sequential manner by application 
of RLS (recursive least square) technique. Because the error matrix G depends on a discrete time k. during a second 
and a subsequent iteration of the decoding operation, it becomes necessary to update the adaptive filter coefficient 
w (k from symbol to symbol, and as mentioned previously, it is preferred to determine the adaptive filter coefficient 
w n (k) by using the channel matrix H. B l 

Fourth aspect of the invention (channel estimation) 

[0069] Using not only known information such as a unique word in the iterative channel estimation, but also using a 
hard decision value of information symbol, in particular one which is likely to be certain as a reference signal is applicable 
not only in the described multiple input multiple output turbo-reception method, but also generally to a turbo-reception 
metttod hn wh.ch a channel (transmission path) of a received signal is estimated from the received signal and the known 
signal, the estimated channel value is used to process the received signal and to decode it. and the decoded signal is 
used in iterating the processing according to the estimated channel value and the decoding processing upon the same 
received signal. 

[0070] Fig. 8 shows an example in which a hard decision value of information symbol is used in a channel estimation 
turbo^quahzer 41. The turbo-equalizer 41 determines a linear equalization filter coefficient in accordanc with an 



16 



EP 1 233 565 A2 



stimated channel value. Th received signal is processed by such a linear equalization filter, th pr cess d signal is 

S?t?2.?SSr? f° dedsi ? nal ' suti,izedinth ite™«v processing of the sam received signal. A received signal 

li, J22£I? UrtX> t. qUa ,Z H r 41 and is also f 0 to a ch annel estimator 42 wh re a channel value (transmission 

SSSTTT^" e l - °" ^ toSiS ° f 0,6 reCeiVed Si9nal and a uni 0 ua word a ™rnory 29. The 
^?„-f f( I ' S 1 J i° 80 eaua,ization Posing in the turbo-equalizer 41 in accordance with the estimated 
channel valuejwd is then subject to a decoding processing, whereby decoded data c'(i) and a soft decision value b' 

ih^h f / ? e Hfl aS,0n ValUe b ' (0 iS inpUt t0 a symbo1 desi 9"ator 43. If its soft decision value^fl) has an 
absolute magnitude wh.ch .s equal to or greater than a threshold, the corresponding hard decision value Wi) is stored 
in a previous symbol memory 32 as one which is likely to be certain (as having a high reliability) to update it During a 
^nt? T? 9 r T CePh °l pr( T ssina totalization Processing) of the same received signal r(k) the channel es- 
deasion value b(i) of information symbol which is stored in the previous symbol memory 32 

c^LI^^T 1 ™ 41 ^ C ° mpriSe the r6Ceiver shown in Fi9 - 1 - for exam P |a - *°m which the iterative 
channel estimator 28 the unique word memory 2 9 and the previous symbol memory 32 are removed. It may comprise 

soluT 9 ' ' 9£ " n " 8 SO ' Uti0n ° f *• eqUati ° n (19) Wi " be in followin 9 form accordi "9 to «he Wiener 



w(k)=E|y(k)y' H (k)]E[b(k)y(k)] 

=[HA(k)H+a 2 I]h (50) 

where H is as defined by the equation (8). and 

h^[H(Q-1) H(0)] T 

where H( ) is as defined by the equation (5), a 2 =E[||v|| 2 ] (variance of noises), and 

A(k)=diag[1-b' 2 (k+Q-1) 1 1-b' 2 (k-Q+1)] 

!2u n f ,SO in *Z Sh ° Wn Fi9 - 31, 0,6 Channel matrix H < > is estimated - and *° channel 

™i T ?J dete . rmine 106 equalization filter coefficient w(k). the received signal is filtered according to the 

I ormtS f " d K t PrOCe ^ d ° UtPUt iS SUbjeCt 8 deCOdin9 Posing. Accordingly, by using a hard decision 
information symbol of a h.gh rel.ab.Hty in the channel estimation during the iterative recessive processing, there can 
be obtained a more correct channel estimation. 

55? t 09 ' 9 Sh *° WS an axample of a torbo-receiver in which the iterative channel estimation method is applied in 

L ' 6 at ° r 42 ' DUnn9 a " initial reception - a channel value is estimated °" ^e basis of the received 
signal r(k) and . unique wort ,n the channel estimator 42, and a compensation for a phase rotation which each symbols 
hasexpenenced on a transmission path as well as RAKE synthesis take place in accordance with the estimated channel 
I« ™ n ♦ l^f T^V ? processor 45 ' or a tim e diversity processing takes place to be delivered to a turtoo-decoder 
to ™^°h ? ? t T deCOded d3ta C ' (i) and 3 S ° n decision value h'W-The soft decision value b-(i) is input 

t 1TE L >• a H d 88 deSCribed eXamPleS> 3 hafd deCfei0n Value b « of a " informati °" ***** whfch 

is likely to be certa.n is stored in a previous symbol memory 32 by updating it. During a second and a subsequent 

£££ -T P r Pr0Ce l S ' ? 9 ° f RAKE reCepti ° n tu *o-decoding, not only the unique word but also the hard decision 
™™ '7°™ " which is obtained during the previous iteration are used in the channel estimator 42 for 

P n( £°, J- estimation, thereby making the channel estimation more accurate to improve the quality. 

ESL m S' ° r a " eX8mple of a torbo-receiver in which the described iterative channel estimation method is 
apphed n the iterative reception us,ng an adaptive array antenna. A received signal r (k) is received by an adaptive 
array antenna receiver 47. and is then branched into a channel estimator 42 where a channel estimation takes place 
in accordance with the receded signal in combination with a unique word. The estimated channel value is used to 

T TZZZTJ^fJ" aPP i6d *° e3Ch a " tenna element 0r spending reception paths in an array weight 
d trm.n.ngun.tsothatth pnnapal beam of the ant nna directivity respons ofth adaptive array antenna r ceiver 
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47 is directed toward th oncoming direction of an intended wave while null is directed toward the ncoming dir ction 
of an interfenng wav , and such weights are applied to applicable locations. A received output from the adaptiv arrav 
antenna receiv r 47is fed to a turbc-decod r 46 tab decoded, whereby decoded data c'(i) and a soft decision value 
b (,) are del* red. Th soft decision value b'(i) is input to a symbol designator 43. and a hard decision value which is 
likely to be certain is stored in a previous symbol memory 32 by updating it. During a second and a subsequent iterative 
reception process.ngof the adaptive array antenna receiver47-turtKMlecoder 46. the hard decision value of information 
symbol which is obtained during the previous iteration is used in the channel estimation in the channel estimator 42 
together with the unique word. In this manner, the channel estimation takes place more correctly, resulting in a more 
accurate control of the antenna directivity response and improving the quality. 

[0075] The turbo-equalizer 41 shown in Fig. 8 is schematically shown in Fig. 11A. and as shown, it comprises a 
series connection of a soft input soft output (SISO) equalizer 41a and an SISO decoder41b. and the operation iterates 
bebveen me equalizer 41a and the decoder 41 b. The turoo-decoder 46 shown in Figs. 9 and 1 0 is schematically shown 
in Fig. 11B. and as shown, it comprises a series connection of a SISO decoder 46a and an SISO decoder 46b and 
an iterative decoding operation takes place between the decoders 46a and 46b. The turbo-decoder 46 shown in Figs 
9 and 10 may comprise a single SISO decoder. ' 
[0076] The examples shown in Figs. 8 to 10 are collectively illustrated in Fig. 12. Thus, a received signal is initially 
processed in an iterative receiver (turbo-receiver) 49 in accordance with a channel value which is estimated by a 
channel estimator 42. the processed signal is decoded, and a result of the decoding operation is provided in the form 
of decoded data (symbol) c-(i) and its soft decision value b'(i). The soft decision value b'(i) is compared against a 
Oireshold value in a symbol designator 43 to determine whether or not a corresponding decoded data cYi) (symbol 
hard decision value) is likely to be certain. If it is determined to be likely to be certain, the hard decision value is stored 
in a previous symbol memory 32 by updating it. and during a second and a subsequent iteration of the processing 
using the estimated channel value and the decoding processing, the symbol hard decision value obtained during the 
previous iteration is used in addition to known information such as unique word in the channel estimation which takes 
place in the channel estimator 42 to provide a more correct channel estimation. 

[0077] Fig .13 shows an exemplary processing procedure of an iterative turbo-reception method which also uses a 
symbol hard decision value. At step 81. a channel value is estimated on the basis of a received signal and a known 
signal. At step 2, an examination is made to see if this is the first one of the iterative processing, and if it is, the channel 
value which is estimated at step S1 is used to process the received signal at step S3, and then the decoding processing 
takes place to determine a symbol hard decision value and a soft decision value. At step S4, a symbol hare decision 
value whjch corresponds to the symbol soft decision value and which is likely to be certain is extracted, and at step 
S5. the extracted symbol hard decision value is used to update a previous symbol hard decision value which is stored 
in the memory 32. At a step S6. an examination is made to see if the number of decoding processings has reached a 
given number, and ,f not. the operation returns to step S1. When it is found at step S2 that the iterative processing is 
not for the first time, a previous symbol hard decision value is read out from the memory 32 at step S7. and it is used 
together with an information symbol of a received signal in performing a channel estimation, whereupon the operation 
ir snsTGrs to stsp S3. 

[007$] Again, as mentioned previously in connection with steps SV to S4' with reference to Fig. 6, the known signal 
may not be used during a second and a subsequent iteration. 

[0079] In the example shown in Fig. 10, RAKE synthesis processor 45 may be inserted between the adaptive array 
antenna receiver 47 and the turbo-decoder 46, as indicated by broken lines. In this instance, a channel estimation 
which is necessary for the compensation of a phase rotation of each symbol and the RAKE synthesis in the RAKE 
synthesis processor 45 may be served by the channel estimator 42, or may take place separately. 

Noise other than white Gaussian noise 

[0080] In the embodiment of the turbo-reception method (according to the first aspect of the invention), the embod- 
iment according to the second aspect of the invention which takes an error correction into consideration and the em- 
bodiment of the turbo-reception method characterized in its channel estimation method (according to the fourth aspect 
of the invention), the processing took place on an assumption that the noise is white Gaussian noise. Specifically v 
M appearing on the right side of the equation (29) indicating a received signal r m (k) from each antenna is assumed 
to be white Gaussian noise. What is meant by white Gaussian noise is a signal which follows the Gaussian distribution 
and have statistical features expressed as follows: 
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10 



B lV Jk >V Jk^f 2 fOT 1 = ° 
0 for q*0 

where E| ] represents an expected value and a variance. White Gaussian noise may be exemplified by thermal 
no.se which is generated in an antenna element. What is influenced by the assumption of white Gaussian noise is & 
^ 1 a P peann 9 eauation O*) whi °" determines the filter coefficient w n (k) or the equation (50) which deter- 

ZZ£2E* Wn(k) " For example ' Wn(k) appearin9 in *• equaUon (34) is calcul * ted thr ° ugh me 
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w n (k)=(HG(k)H H +E[n(k)n H (k)])-^ 
=(HG(k)H H +0 2 I)' l h 

where v m (k) is calculated as E[n (k)nH(k)) = a*! by the assumption of the white Gaussian noise having a variance o* 
The channel matnxH which is estimated by the iterative channel estimator 28 (Fig. 1)or42(Fig 12) <y2 and the error 

sssizsi izssszzz* basis of the a priori ,o9 - nke,ihood vaiue are substituted int ° * ™ to 

r^ 1 i a rth^ e a n ° iSe « V ^ (k f iS ^ ite GaUSSia " n0iSe> E[n (k) nH(k)1 =o2 ' does not a PP ,v - Accordingly, in order 
to calculate the filter coefficient w n (k), it is necessary to estimate an expected value (covariance) matrix E[n (k) nH(k)] 

InriTT' V f sa P a ^ted method^uch method will now be described. A covariance matrix for the noise 

hST, ! ^ * US E[ " <k) n (k)L yW =H - B ( k ) + "< k > ^ equation (22) is modified into n(k)=y(k) 
- H-B(k) and is substituted into the covariance matrix U, as indicated below. 

U=E[n(k)-n H (k)] 
=E[(y(k)-H-B(k))-(y(k>HB(k)) H ] 

[0082] If we can assume that a matrix y(k) is available from a received signal, an estimate H of a channel matrix H 
from the channel estimate and B(k) is available from a reference signal, it is possible to estimate the matrix U according 
to the time average method as follows: 

Tr (y(k)-H B(k)) (y(k)-H B(k)) H (51) 
where Tr represents the number of reference symbols. 

[0083] During an iterative channel estimation which takes place in the iterative channel estimator 28 or 42 the co- 
vanance matnx 0 is estimated using the channel matrix H together with the equation (51). A procedure therefor is 
illustrated in Fig. 14. Fig. 14A shows a unique word and information symbol series which occurs in one frame of a 
received signal, and Fig. 14B shows the initial processing as well as a subsequent processing. During the initial process- 
ing on y the unique word is used as a reference signal, initially estimating the channel matrix H. Subsequently U is 
estimated using the unique word and the channel matrix estimate H according to the equation (51). Using the estimate 
U and H, the filter coefficient w„(k) is calculated as follows: 8 



55 



W„(k)=(HG(k)H H +0)- 1 h (52) 

fnfonStion sjmbof ""^ " ' 0 **** 3 ** equa,ization upon 0,6 received si 9" a| . thus estimating transmitted 
[0084] During a second iteration, the unique word as well as one of information signals (*) estimat d during th initial 
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qualizabon which is d termined to be life lyt b certain according to the threshold value are both used as reference 
signals to stimate H again with the same procedure as;us d during the initial processing, thus stimating U again 
As this operation is repeated, the channel matrix stimate H becomes more accural withth it ration, and th stimat 
of U becom s more accurate, thus improving th accuracy of the filter coefficient w n (k) to improve the response of the 
equalizer. 

[0085] In this manner, a turbo-reception when a noise other than white Gaussian noise is contained in a received 
signal is made possible. 

[0086] A functional arrangement in which a linear equalization is performed by estimating a covariance matrix U of 
a noise contained in a received signal is shown in Fig. 15 as an example in which a log-likelihood ratio A« [b«(k)] is 
obtained as an equalization output from a multiple output equalizer 31 shown in Fig. 2 for a received signal correspond- 
ing to a transmitted signal from a first transmitter. Parts corresponding to those shown in Fig. 2 are designated by like 
reference numerals in Fig. 15 as used in Fig. 2. 

[0087] A unique word from a unique word memory 29 or a previous symbol hard decision from a previous symbol 
memory 32 which is likely to be certain is input to a reference matrix generator 31 9, which then generates a refe/ence 
matnx B(k) according to the equations (25) and (26). The reference matrix B(k), an estimation channel matrix H from 
a channel estimator 28, and a received matrix y (k) from a received matrix generator 31 1 are supplied to a covariance 
matrix estimator 321 , which then calculates the equation (51 ) to obtain an estimated matrix 0 for a covariance matrix U 
[0088] Soft decision transmitted symbols b^k) ,-.,b n '(k) from a soft decision symbol generator 313-1 are input to an 
error matrix generator 322-1 , where an error matrix G,(k) corresponding to the square error of the channel estimation 
is generated according to the equations (35), (36) and (37). The error matrix G^k), the estimated covariance matrix 
U and the estimated channel matrix H are supplied to a filter estimator 323-1, where the equation (52) is calculated 
to estimate a filter coefficient Wl (k). The filter coefficient w,(k) and the difference matrix y'(k) from a difference calculator 
316- are fed to an adaptive filter 318-1 where a filter processing w^kjHyfk) is applied to y'(k), and its result is delivered 
as a log-likelihood ratio A 1 [b^k)]. 

[0089] When reflecting an error correction recording result into a signal being detected, a function calculator 331-1 
as shown in Fig. 7A is provided, as indicated by broken lines in Fig. 15, to calculate f(b' n (k)), and an interference replica 
matnx generator 314-1 uses the equation (43) rather than the equation (31), and a error matrix generator 322-1 uses 
the equation (48) rather than the equation (37). 

[0090] The procedure shown in Fig. 14B is shown as a flow chart in Fig. 16. Specifically, at step S1 , a received signal 
r(k) and a known signal (for example, a unique word) are used to estimate a channel matrix H f and at step S2 an 
examination is made to see if this processing is for the first time, and if it is, the known signal, the estimated channel 
matnx H and the received signal r(k) are used to calculate the equation (51) to determine an estimated covariance 
matnx 0 at step S3. A 

[0091] At step S4, the estimated channel matrix H, the estimated covariance matrix 0 and an error matrix G(k) which 
compnses symbol soft decision values are used to calculate the equation (52) to estimate a filter coefficient w n (k) 
[0092] At step S5, the estimated channel matrix H and the filter coefficient w n (k) are used to equalize the received 
signal or to calculate the equation (27) to determine w„h (k) y'(k) to obtain a log-likelihood ratio A-j [b n (k)], subsequently 
performing a decoding process to estimate a hard decision value and a soft decision value of a transmitted symbol 
[0093] Purpose of step S6 is to determine a symbol hard decision value which corresponds to a symbol soft decision 
value which exceeds a threshold value and which is likely to be certain (or having a high reliability). This symbol hard 
decision value is used to update a symbol hard decision value which is stored in a previous symbol memory 32 Sub- 
sequently, an examination is made at step S8 to see if the number of the decoding processings has reached a given 
value, and if not, the operation returns to step S1. However, if a given number is reached, the processing upon this 
received frame is completed. 

[0094] If it is found at step S2 that the iterative processing is not for the first time, namely, for a second or a subsequent 
iteration, a symbol hard decision value is read from the previous symbol memory 32 at step S9, and is used together 
with information symbol in the received signal to estimate the channel matrix H, subsequently transferring to step S3 
[0095] Again, by changing the steps S1 and S2 in the similar manner as steps SV to S4' shown in broken lines in 
Fig. 6, it is possible to avoid the use of the known signal during a second and a subsequent iteration. When it is desired 
to reflect an error correction recording result into a signal being detected, a function f(b' n (k)) may be calculated at step 
S10, as shown in broken lines in Fig. 16, and it may be used to obtain the error matrix G(k). In either instance the hard 
decision transmit symbol may not be used in the estimation of the covariance matrix 0. The capability to estimate a 
covanance matrix U of a noise contained in a received signal which is other than white Gaussian noise is applicable 
to a variety of useful purposes, as will be described below. 

(1) A reception method for a multiple series transmitted signal containing an unknown interfering signal is cited 
As shown in Fig. 30A, it is assumed that an interfering signal i(k) which is unknown to a turbo-receiver such as 
signals from other cells or zones in a mobil communication, for example, to be r ceived by the turbo-receiver in 
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addition to Ns riestransmitt d signal as may b transmitt d from N user transmitters. In this instance, the quation 
(20) can be writt n into th following form: 



r m (k)=I q= o °" 1 2n=i N h mn (q) b n (k-q)+i(k)+v m (k) (20)' 
In this model, putting i(k)+v m (k)=v' m (k), we have 



'mOO^o ^ In=1 N h mn (q)b n (k-q) + v« m (k) (20)" 

Treating v' m (k) as a noise signal other than white Gaussian noise, H and U are estimated in a manner mentioned 
previously, and w n (k) is estimated, and a turbo-reception can be made by iterating an equalization of a received 
signal and an estimation of transmitted symbols. 

(2) In a communication system which employs a transmission/reception separation filter, when an oversampling 
of a received signal is made at a higher rate than a symbol period, there occurs a correlation between noise 
components which are contained in received signals which are sampled at individual times, and this prevents the 
noises in the received signals from being treated as white Gaussian noise. In other words, the equation (20) does 
not apply. Accordingly, an assumption represented as 



E[n(k) n H (k)]=a 2 I 

does not hold. A processing upon a received signal which is separated by the transmission/reception separation 
filter may utilize equation (51) to determine a covariance matrix U, thereby allowing the received signal to be 
properly processed. 

(3) In the described turbo-reception method, the arrangement is such that every multipath component from each 
transmitter (user) on Q paths are synthesized. However, in the event there exists a prolonged delay wave on 
channels (for example, assume that paths include one symbol delay, two symbol delay and three symbol delay 
path and there exist separately a thirty symbol delay path: in this instance, the thirty symbol delay path component 
is treat as a prolonged delay wave), it is possible to prevent the prolonged delay wave from being synthesized 
but to treat it as an unknown interference which can be eliminated by an adaptive filter. When the prolonged delay 
wave component is treated as the interfering signal i(k) in the example according to the first aspect of the invention 
(1), it may be eliminated. 

[0096] In the processing of a received signal containing a noise other than white Gaussian noise, the estimation of 
the covanance matrix U is applicable to a single user turbo-reception method by allowing it to serve in place of o 2 1 in 
the equation (50). In a similar manner, it may be used in a RAKE synthesis processing reception illustrated in Fig 9 
or a turbo-reception using an adaptive array antenna reception shown in Fig. 1 0, irrespective of a single user or multi- 
user application, or more generally, in a channel estimation in a channel estimator 42 and an estimation of the covar- 
iance matrix U during the iterative decoding operation shown in Fig. 1 2. For RAKE reception, only the channel estimation 
may be used. 

Third aspect of the invention (multistage equalization) 

[0097] In the forgoing description, received signals r 1f -,r M are equalized in a multiple output equalizer 31 to deter- 
mine log-likelihood ratios A,[b(k)], A^bfk)], but in a modification (2) of the first aspect of the present invention there 
are provided a plurality of equalizer stages in cascade connection in a manner such that the number of outputs is 
reduced toward a later stage equalizer. By way of example, Fig. 17 shows that the equalizer is divided in two parts 
where a prestage equalizer (multi-user equalizer) 71 cancels an interfering component which is located outside an 
equalization range of a poststage, single user equalizer 21 *. At this end, a preprocessing including a soft interference 
cancel and a linear filtering according to MMSE(minimum mean square error) criteria, for example, is made, and sub- 
sequently, the poststage equalizer 21 * performs an equalization of a single user having a number of paths equal to Q 
[0098] Even when an equalization takes place in a cascade manner and a linear filter is used in a prestage processing 
it is possible to prevent the quantity of calculation from increasing prohibitively. 

[0099] An arrangement of a multiple output turbo-receiver according to an embodiment which is based on this fun- 
damental concept of the first aspect of the invention (2) of the turbo-reception method and an exemplary arrangement 
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of MIMO system which incorporates th pres nt invention is shown in Fig. 18. It is to be understood that parts corre- 
sponding to thos shown in Fig. 1 are d signal d by lik ref rence characters as us d before, without repeating their 
description. (Th same holds true in subs qu nt d scription.) 

[0100] Transmitt d signals from each transmitt r is received through transmission paths (channels) by a turbo-re- 
ceiver 30. The received signal r(k) is input to a multi-user equalizer 71. from which signals from N transmitters are 
delivered in the form of signals Ul (k). u N (k), each of which is provided in the form it is free from interferences by 
signals from other transmitters, and channel values o^k). o N (k) to be input to single user equalizers 21-1 - 21-N 
These SISO equalizers 21-1, -.21-N deliver log-likelihood ratios A, [b,(k)]. -. A, [b N (k)J. Subsequent passing 
remains similar to Fig. 1. However, the channel values o,(k), -. a^k) which are used in the single user equalizers 
21-1. -. 21-N are channel values which are obtained subsequent to the multi-user equalization, and are distinct from 
ttie channel matrix H. Accordingly. ai (k). -. « N (k) are referred to as post-equalization channel information. 
[0101] The operation will now be described more specifically. 

[0102] Equations (23) to (26) are defined in a similar manner as described above in connection with Fig. 1 in con- 
sideration of the number of multipaths (channels) Q. 

[0103] The purpose of the poststage equalizers 21-1 , -.21-N shown in Fig. 18 is to equalize intersymbol interference 
channels by signal symbols of respective users themselves [b n (k). b n (k-1). -. b n (K-Q+1)] (where n=1. -. N). At this 
end. the prestage equalizer 71 operates to eliminate interferences other than [b n (k). b„(k-1), - b n (K-Q+1)l (where 
n=1. N) within y(k). This will be discussed below more quantitatively. 

[0104] Initially, using the a priori information [b n (k)] (where n=1. — ,N) of the equalizer 71 which is fed back from 
ttie decoders 24-1 . -. 24-N. a soft decision transmitted symbol estimate b'(k) is determined according to the equation 
(15). 

[01 05] This soft decision transmitted symbol b' n (k) and a channel matrix H are used to generate a replica H-B'(k) of 
an interfering signal, which is then subtracted from the received matrix y(k). 

y'n(k)^y(k)-HB'(k) (27 )' 
=H(B(k)-B'(k))+n(k) (28 )' 



where 



B'(k)=[b' T (k+Q-1)...b' T (k)...b ,T (k-Q+1)] T (29)' 
b'(k+q)=[b* 1 (k+q)b' 2 (k+q)...b' n (k+q)...b' N (k+q)] T :q=Q-1 1 (53) 



b , (k+q)=[b , 1 (k+q)b' 2 (k+q)...0...b' N (k+q)] T :q=0 -Q+1 



(54) 



it being noted that b'(k+q) has a zero element at n-th position. 

[01 06] The operation of subtracting the interference in this manner will be referred to hereafter as a soft interference 
cancel. Assuming that a replica of the interfering signal is generated in an ideal manner, it will be seen that y' (k) 
resulting from the subtraction can have only the symbol b n (k) of the n-th user and an intersymbol interference compo- 
nent caused by the symbol [b n (k-1), b n (k-Q + 1)] of the n-th user itself which results from putting the n-th element of 
b (k+q) equal to 0 for q=1, ■« -Q+1 in the equation (54). 

[0107] In effect, a contribution from the signal of the n-th user (transmitter) to the received matrix r(k) is only that 
resulting that from the symbol [b n (k), b n (k-1), b n (k-Q+1)]. However, it will be understood from the definition of the 
receded matrix y(k) given by the equation (21) that a contribution from the signal of the n-th user (transmitter) within 
the received matrix y(k) which results as synthesis of the multipaths will contain, when based on the k-th symbol b 
(k), intersymbol interference components caused by future symbols [b n (k+Q-1), b n (k+Q-2), ■» b n (k+1)] Thus the in- 
terference replica includes interference components from the future symbols. In this respect, the difference matrix y 
(k) defined by the equation (27)* is distinct from the difference matrix y'(k) defined by the equation (27) 
[0108] Accordingly, a next step in the prestage processing in the equalizer 71 is to eliminate the residue interference 
which remains after the soft interference cancel, namely, a residual interference component which results from an 
imperfect synthesis of the interference replica H-B'(k) and interference compon nts b twe n future symbols from / 
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(k) using MMSE (minimum mean square rror) crit ria lin ar filter. In oth r words, this elimination takes place by an 
arrang m nt such that a filtering of y' n (k) using th filter characteristic w n (k) as indicat d by the equation (55) is qual 
to a sum f the symbols [b n (k), b n (k-1), b n (K-Q+1)J ach multiplied by channel values a 1n (k), a^k), • -.a^k). 

W n H (k)y n (k) ~ 2q=o Otq (k>b n (k-q) 

=a n H (k>b n (k) (55) 



Accordingly, what is required is to calculate the equation (55) by determining the filter characteristic w n (k) and the post- 
equalization channel value (channel information) o^k). The derivation of w n (k) and a^k) will be described. It is to be 
noted that while the filter characteristic w n (k) is distinct from the filter coefficient w n (k) given by the equations (32) and 
(34), similar denotations will be used for purpose of convenience. 
[0109] Desired solutions are defined as solutions of the following optimization problem: 



(w n (k),a n (k))=arg min||w n H (k).y' n (k)-a n H (k) b n (k)f (5 6) 

provided a 1n (k)=1. 

[0110] In other words, w n (k) and a n (k) which minimizes the right side of the equation (56) are determined. The 
constraint requirement a 1n (k)=1 is added in order to avoid solutions which result in a n (k)=0 and w n (k)=0. While solutions 
can be obtained under the constraint requirement, 



IK(k)|| 2 =1 

a solution for a 1n (k)=1 will be described below herein. For brevity, the problem will be redefined. Namely, the right side 
of the equation (56) is defined as m n (k) which is minimized in terms of w and a. 



m n (k)=arg min||m n H (k)-z n (k)f (57) 
35 provided m n H (k) e MQ+1 = -1 (which is equivalent to a 1n (k) =1) and where 

m n (k)-[w n T (k),-a n (k) T ] T (58) 

40 

z n (k)-[y n T (k), b(k) n T ] T (59) 



45 



e MQ+1 =[0...1...0] T 



(60) 



it being understood that e MQ+1 has "1 n element at (MQ+1-th position). 

[0111] A solution of the optimization problem is given as follows according to Lagrange's method of indeterminate 
coefficients described in literature [2], S. Haykin, Adaptive Filter Theory, Prentice Hall, pp. 220-227; 

50 

mnW^Rzz^eMQ+i^MQ+^Rzz'^eMQ+i) (61) 



where 

55 



R zz =E[z n (k)z n H (k)] ( 6 2) 
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E[A] representing an xp cted value (a mean value) of A. 



= E 



H-A n (k)-H H +<x 2 l H* 



I 



(63) 



A n (k)=diag[D n (k+Q-1),...,D n (k) D n (k-Q+1)] 

where I represents a unit matrix, and o 2 noise power (a covariance of white Gaussian noise). 



(64) 



hn(Q-l) 0 0 
h n(Q-2) hn,(Q-l) 0 



Vo) 



Ml) 



0 
0 



(65) 



D n (k+q)=diag[1-b* 2 1 (k+q) 1-b' \ (k+q) 

1-b ,2 N k+q]] :q=Q+1 1 



(66) 



D n (k+q)=diag[1-b' \ (k+q) 1 1-b^ (k+q)] 

:q=0 -Q+1 



(67) 



where diag represents a diagonal matrix (all elements other than those located along the diagonal of the matrix being 
zeros). Thus, if the channel matrix H and a 2 are known, m n (k) can be determined according to the equation (61). 
Accordingly, w n (k) and ot^k) can then be determined according to the equation (58). 

[0112] Using the filter characteristic w n (k) which is calculated in this manner, y' n (k) is filtered according to the following 
equation; 



"n(k)=w n M (k) y' n (k) 



(68) 



where H represents a conjugate transposed matrix. 

[0113] These filtered n processed results are fed to corresponding equalizers 21-n which follow. In this manner, a 
received signal u n (k) which corresponds to the left side of the equation (1) from the n-th user is obtained, a^k) which 
corresponds to the channel value h mn (q) on the right side of the equation (1) is obtained, and the equation (55) which 
corresponds to the equation (1) is determined. Accordingly, ccjk) is applied as an equalizer parameter (channel value) 
to a succeeding equalizer 21 -n. This completes the prestage processing by the equalizer 71 . 
[0114] A processing which takes place in the succeeding equalizer 21-n and thereafter will now be described. As 
mentioned previously, because the equation (55) corresponds to the equation (1), the operation which takes place in 
the equalizer 21-n for every user may proceed in the similar manner as that of the equalizer 21 shown in Fig. 31 , and 
will not be repeated again since it is disclosed in literature [1]. Each equalizer 21-n receives u n (k), ctn(k) and a priori 
information 7^ [b n (k)] from the decoder 24-n as defined above, and calculates and delivers outputs a log-likelihood 
ratio A«j ( LLR) representing the ratio of the probability that each encoded bit is +1 and the probability that it is -1 , as 
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follows: 



A ru,^-,^ Pf P>n(k) = ^KOO. k ■ 0, -,B] 



Prtb n {k) = -1|u n (k).k = 0.-,B] < 69 > 



-Xi[b n (k)] + X|tb n (k)] (70) 



where ^ [b n (k)] represents an extrinsic information fed to a succeeding decoder 24-n, and 7^ [b n (k)] a priori information 
applied to the equalizer 31 . Decoder 24-n calculates the log-likelihood ratio A 2 according to the following equation: 

a rh /:u-«~- P ^i) = +1l*iP>n (i)Li = 0/",B] 

A2lbM ~ l09 PriW (71) 



-Mb n (i)] + X?[b n (i)J (72) 

where ^ [b n (i)] represents an extrinsic information applied to the equalizer 71 and the equalizer 21 during the iteration 
and the [b n (k)] a priori information applied to the decoder 24-n. The arrangement shown in Fig. 18 performs an 
iterated equalization and decoding to improve an error rate. 

[0115] A functional arrangement of the multi-user equalizer 71 will be briefly described with reference to Fig. 19. A 
received signal from each antenna is processed in a receiver 70 as a matrix rfkHr^k) • r M (k)] f which is used in a 
received matrix generator 311 to generate a received matrix y(k) according to the equation (21) which takes the indi- 
vidual multi-paths (channels) into consideration. 

[0116] On the other hand, the received signal r(k) from the receiver 70, and a known series signal such as a unique 
word series used for channel estimation which corresponds to each transmitter and which is fed from a unique word 
memory 29 are input to a channel estimator 28 in order to estimate a channel matrix H. 

[0117] A priori information ^ p [b n (i)], \f [b N (i)] are subtracted from the log-likelihood ratios A 2 [b^i)], -,A 2 [b N 
(i)] delivered from the respective decoders 24-1, • 24-N to derive extrinsic information [bi(k)], •-, ^ P> N (k)] f which 
is then input to soft decision symbol estimators 31 3-1 , 313-N where soft decision transmitted symbols b'^k), b* N 
(k) are calculated according to the equation (15) and are then input to an interference matrix generator 72. In the 
interference matrix generator 72, a matrix B'(k) of symbol estimates which can be interference signals from other 
transmitters are generated for each n according to the equations (29)*, (53) and (54). A product of these N matrixes B* 
(k) and the channel matrix H is generated by other interfering signal estimators 73-1, 73-N, respectively, thus de- 
termining the replica of interfering components H-B(k). 

[0118] These N interfering component replicas H-B(k) are subtracted from the received matrix y(k) in subtractors 
74-1, 74-N, respectively, thus providing difference matrixes y^k), y* N (k). 

[0119] The soft decision transmitted symbols b^k), -.b^fk) are input to an error matrix generator 75 where error 
matrices (k), A N (k) are generated according to the equations (64), (66) and (67). These error matrices, the channel 
matrix H and the noise power a 2 are input to a filter characteristic estimator 76 where the filter characteristic w n (k) and 
the post-equalization channel information a„ are estimated according to the equations (58), (60), (61), (63) and (65). 
These filter characteristics w 1t w N and difference matrixes y'^k), ■•, y' N (k) are multiplied together in filter processors 
77-1, 77-N, respectively, or the difference matrixes are filtered, thus determining a component u^k), u N (k) of the 
received signal for the symbol [b n (k), b n (k-1), b n (K-Q+1)] from each user and for each path, from which interferences 
from other user signals are eliminated. These components are fed, together with the post-equalization channel infor- 
mation cc^k), -, a^k) which are determined in the filter characteristic estimator 76, to the single user equalizers 21 -1 , 
— , 21 -N shown in Fig. 18. 

[0120] A processing procedure for the turbo-reception method according to the first aspect of the invention (2) is 
shown in Fig. 20 where steps corresponding to those shown in the procedure of Fig. 3 are designated by like step 
numbers as used before. However, the calculation of the interference replica matrix B' n (k) which takes place at step 
S4 follows the equations (29)\ (53) and (54). The step S1 3 uses the soft decision transmitted symbol b' n (k) to generate 
the error matrix AJk) according to the equations (64), (66) and (67). Step S14 uses the channel matrix H, the noise 
power a 2 and the error matrix Ap(k) to determine the residual interference eliminating filter w n (k) and the channel 
information o„ according to the equations (58), (60), (61), (63) and (65). Step S15 filters the difference matrix y' n (k) 
according to the residual interference eliminating filter characteristic w n (k) to determine u n (k). At step S16, a single 
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user equalization is applied to each filt red result u n (k) to det rmine the log-likeliho d ratio AJu^k)], which is then 
decoded at st p S10. In oth r respects, the procedure is similar to that shown in Fig. 3. 

[0121] Inth f rgoing description, the extent of qualization in the poststage equalizer 21 -n is defined as a zone for 
mt rsymbol interference by the symbol [b n (k), b n (k-1). b n (K-Q+1 )] (where n=1 , N), but such extent of equalization 
is adjustable. For example, when Q has a very high value, a computational load on the poststage equalizer 21 -n will 
be much greater. In this instance, the extent of equalization by the poststage equalizer 21 -n is chosen to be Q'<Q while 
the prestage equalizer 71 may be rearranged so as to eleminate intersymbol interferences between the signal of the 
same user outside the zone [b n (k). b n (k-1), ■-, b n (K-Q' + 1)] (where Q'<Q and n=1, N). Such modification will be 
descnbed later. When the equalization is divided into the prestage and poststage. a previous symbol memory 32 may 
be provided, as indicated in broken lines in Fig. 19, in association with the channel estimator 28 so that the hard decision 
transmitted symbol b n (k) may also be used in the estimation of the channel value, thus permitting the accuracy of 
estimation to be improved. 

[0122] In the example shown in Fig. 17, the prestage multiple output equalizer 71 equalizes transmitted signals in 
N senes so that interferences from other series may be separated to provide signals u n of N series and the post- 
equalization channel information and subsequently, the signal u n of each of N series is processed by the post stage 
single user equalizer 22-n so as to eliminate intersymbol interference of the same transmitted signal. In this manner 
the equalization takes place in two stages which are in cascade connection. However, a cascade connection of three 
or more stages may also be used. 

[0123] By way of example, Fig. 21 shows that a received signal r m of M series with respect to a transmitted signal 
of N senes is input to a first stage equalizer 81 to provide an equalized signal series er^k) for 1st to U-th transmitted 
senes from which interferencs by (U+1)th to N-th transmitted series is eliminated and its associated post-equalization 
channel information ea(k), and an equalized signal series er 2 (k) for (U+1)-th to Nth transmitted series from which 
interferences by 1st to U-th transmitted series are eliminated and its associated post -equalization channel information 
eo2(k) while a second stage includes equalizers 82-1 and 82-2. er^k) and ect^k) are input to the equalizer 82-1 where 
they are equalized to provide an equalized signal series er 3 (k) for 1st -U r th transmitted series among the 1 st to Uth 
transmitted series from which interferences by (U^IJ-th to U-th transmitted signal are eliminated and its associated 
post-equalization channel information ect 3 , an equalized signal series er 4 (k) for (U^^-th to U 2 -th transmitted series 
among the 1st to U-th transmitted signal from which interferences by the 1-st to U r th transmitted series and U 2 -th to 
U-th transmitted signal are eliminated and its associated post-equalization channel information ea 4 (k), and an equalized 
signal senes er 5 (k) for (U 2 +1)-th to U-th transmitted series among 1-st to U-th transmitted series from which interfer- 
ences by 1-st to U 2 -th transmitted series are eliminated and its associated post-equalization channel information ea 5 (k) 
[0124] Similarly, the equalized signal series er 2 (k) and the channel information ea^k) are input to the equalizer 82-2 
in the second stage to provide an equalized signal series er 6 (k) and an associated post-equalization channel information 
ea^k), and an equalized signal series er 7 (k) and its associated post-equalization channel information ea 7 (k): When 
N=5, equalizers 83-1 to 83-5 in a third stage represent single user equalizers shown in Fig. 18. Alternatively, an input 
equalized signal to the equalizer 83-3 may comprise two transmitted signals, and the equalizer 83-3 can eliminate 
mutual interferences between the two transmitted signals to provide a set of two equalized signals and their associated 
post-equalization channel information, which are in turn supplied to and equalized by next succeeding single user 
equalizers 84-1 and 84-2. As a further alternative, the equalizer 83-4 may receive the equalized signal er 6 (k) and the 
channel information eaefk) to eliminate mutual interferences for each of the transmitted signals which constituted the 
equalized signal er 6 (k), such as each of three transmitted signals, each of which may be interfered by two other trans- 
mitted signals as well as an intersymbol interference of each transmitted signal due to multipaths. As another alternative 
one or more of the equalizers 82-1 and 82-2 in the second stage may be arranged so that an equalized signal may be 
simultaneously obtained for each of a plurality of transmitted signals. 

[01 25] What has been described above can be generalized by stating that equalizers in a first stage deliver a plurality 
of equalized signal series and a set of post-equalization channel information, and one or more equalizers may be 
provided in each of one or a plurality of stages which are in cascading connection for each equalized signal series and 
its associated set of post-equalization channel information so that an equalized output or a log-likelihood ratio A« [b 
(k)] may be delivered finally for each of the 1 st to N-th transmitted series. " 
[0126] When the equalization takes place in multiple stages which are in cascade connection, it is preferred that the 
number of paths Q for which interferences are to be eliminated be reduced toward a later stage so that the quantity of 
calculation can be reduced. In this instance, an interferring component from a path which disappears in a later stage 
be eliminated in an immediately preceding equalizer stage. 

[0127] An equalization processing which occurs in the arrangement of Fig. 21 when a first stage equalizer 81 deals 
with N transmitted signals, each of which has a number of multi-paths equal to Q to provide a group of equalized signal 
senes er^k) comprising U transmitted signals and their associated post-equalization channel information ea^k) and 
a later stage equalizer 82-1 performs an equalization for a number of multipaths equal to Q' for each transmitted series 
will be d scribed. 
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[0128] In a similar manner as d scrib d above in conn ction with the embodiment shown in Figs. 18 and 19. an 
interference matrix generator 72 g n rat sanint rference matrix B'(k), but th quationsus d are changed from the 
equations (53) and (54) to quations (53), (54)' and (73) as follows: 

b t (k+q)=[b' 1 (k+q)b' 2 (k+q)...b' n (k+q)...b' N (k+q)] T :q=Q-1,...1 ( 5 3) 
b'(k+q)=[0...0 b' u+1 (k+q)...b' N (k+q)] T :q=0 -Q'+1 (54)' 

b'(k+q)=[b* 1 (k+q)b' 2 (k+q)...b' n (k+q) 

-b' N (k+q)] T :q=Q' -Q+1 (73) 

The equation (54)' is intended to provide symbols for a first to U-th transmitted series themselves, and to provide an 
equalization except for intersymbol interference caused upon each series by itself and relative to each other due to 
multipaths Q' while the equation (73) is intended to eliminate the intersymbol interferences upon the 1st to U-th trans- 
mitted series themselves and relative to each other due to (Q'+1 )-th to Q-th path in as much as the number of multipaths 
is reduced to Q' in the poststage equalization. 

[0129] The interferences matrix B'(k) which is obtained in this manner is used to generate an interference signal 
replica H-B'(k) which is then subtracted from a received matrix y(k) as follows: 

Y'g (k)^y(k)-H B'(k) ( 27)" 

=H(B(k)-B*(k))+n(k) (28)" 

[01 30] This operation of subtracting the interference is referred to hereafter as a soft interference cancel. Assuming 
that a replica H-B'(k) of an interfering signal is generated in an ideal manner, it will be seen that y'(k) can only have 
signal components for symbols [b n (k), b n (k-1), b n (k-Q'+1)] (where n=1 to U) for 1st to U-th transmitted symbols. 
[0131] A residue of interferences which remain after the soft interference cancel is eliminated with a linear filter of a 
MMSE criteria in a similar manner as mentioned previously. In this instance, the equation (55) is replaced by an equation 
(55)' indicated below. 

w g H (k)-y g (k)«I n=1 u I q=0 Q '" 1 a nq (k) b n (k-q)=a g H (k) b g (k) (55)' 

where 

a g( k ) = KoO<) <*1.Q'-1 (k) a uo (k) a UtCM (k)] T (55-1) 

b g (k)=[b 1 (k) b^k-Q'+l) bu(k) b^k-Cr-M)] 1 (55-2) 

[01 32] The derivation of w g (k) and ctg(k) takes place in a similar manner as described previously to determine w fl (k) 
and ctg(k) which minimizes the right side of the following equation, which stands for the equation (56): 

(w g (k),a g (k))=arg min||w g H (k) y' g (k)-a g H (k) b g (k)|| 2 (5 6)' 

provided a 10 (k)=1. 

[0133] The constraint requirement is added in order to avoid solutions which may result in ctg(k)=0 and w g (k)=0. 
While a constraint requirement thatUog(k) |f=1 may also be used to provide a solution, in the description to follow the 
probl misr writt n as follows for 0 (k)-1 : 
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m g (k)=arg min||m 9 H (k) z g (k)|| 2 



(57)' 



provided m g H (k)e MQ . +1 = -1 



m g (k) s [w g T (k).-o g T (k)J T 



*oW-&g T W*W B T | T 
e M Q' + i=[0...1...0] T 



(58) ' 

(59) " 

(60) ' 



it being understood that e MQ . +1 has "1" element at (MQ'+i)-th position. 

[0134] A solution of the optimization problem is given as follows according to Lagrange's method of indeterminate 
coefficients disclosed in literature [2]; 



m g (k)=-R zz " 1 •e MQ . +1 /(e MQ . +1 H -R 2Z - 1 e MQ , +1 ) 



R 22 = E[Z g (k)zy(k)] 



(61)' 



(62)' 



= E 



HA n (k).H H +a 2 l Hf 



H g I 



(63)' 



A n (k)=diag[D„(k+Q-1) D n (k) D n (k-Q+1)] 



(64') 



h l(Q-l) 


0 0 




0 0 




h l(Q-2) 


0 


•• h u( Q_ 2 ) 


0 




hl(Q-3) 


: hl(Q-l) • 

• • 


- h U( Q_3) 

• • j 


hlJ(Q-l) 

■ • 
• • 


(65)' 


hl (0 ) 


— hi(Q'_D • 


h U( o) 


— hu(Q'-l) 





D n (k+q)=diag[1-b ,2 1 (k+q) l-b^k+q),. 

l-b^fk+q)] :q=Q+1 1 



D n (k+q)=diag[1 1,1 -b^ (k+q), .... 



(66) 



1-b^ N (k+q)] :q=0 -Q'+1 



(67)' 
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D n (k+q)=diag[1-b' 2 1 (k+q) l-b'Vk+q) 

1-b^ (k+q)] :q=Q\...-Q+1 ( 74) 

Thus when channel parameters are known, m g (k) can be determined according to the equation (61)\ In addition, w g 
(k) and Og(k) (= ecc^k)) can be determined according to the equation (58)'. Such calculations can be made in the filter 
characteristic estimator (76) shown in Fig. 1 9, for example, and the following equation is calculated in the filter processor 
77-1 for purpose of filtering: 

er, (k)=w g H (k)y g (k) 

This equalized output en, (k) and the post-equalization channel information ea 1 (k)=a a (k) are fed to the poststage equal- 
izer 82-1. 

[0135] When there are five transmitted series (users), for example, which are divided into a group of three transmitted 
series (users) and a group of two transmitted series (users), in the manner mentioned above, the above algorithm is 
carried out with U= 3 and 2, and the two equalized outputs er^k), ea^k) and er 2 (k), eo^k) are input to equalizers 
which are designed to deal with the three transmitted series and the two transmitted series, respectively, thus obtaining 
an equalized output for each transmitted series. 

[0136] Reflecting an error correction decoding result for a signal being detected into a soft decision transmitted 
symbol in the manner mentioned above is also applicable to a single user turbo-equalizer receiver shown in Fig. 8, 
RAKE synthesis turbo-receiver shown in Fig. 9, a turbo-receiver including an adaptive array antenna receiver shown 
in Fig. 10, a generalized turbo-receiver including a channel estimator 42 shown in Fig. 12. 

[0137] In Figs. 13, 14 and 15, a symbol hard decision value which has been determined as being as likely to be 
certain is also used as a reference signal in the estimation of the channel matrix H and the covariance matrix 0 during 
a second and a subsequent iteration. However, during a second and a subsequent iteration, only the unique word may 
be used as a reference signal to utilize the equation (51) to estimate the covariance matrix 0 while omitting the esti- 
mation of the channel estimation and the estimation of the covariance matrix 0 which utilize the symbol hard decision 
value. 

First aspect of the invention (2) (parallel transmission) 

[01 38] There is a proposal that information series c(i) from a single user be transmitted in a plurality of parallel series 
in order to achieve a high rate transmission with a high frequency utilization efficiency. An embodiment of a turbo- 
receiver incorporating the present invention which may be used for such transmitted signal will now be described. 
[0139] Referring to Fig. 22, where parts corresponding to those shown in Fig. 1 are designated by like reference 
characters as used before, on the transmitting side, a modulated output signal bfl) from a modulator 13 is input to a 
series-parallel converter 14 where each symbol b(j) is sequentially distributed into N series. It is assumed that there 
are series signal b^k), b N (k), the number N of which is equal to or greater than 2. While not shown, these signals 
are transmitted from N antennas after convention into radio frequency signals. 

[0140] These N series signals are propagated through channels (transmission paths) to be received by the turbo- 
receiver according to the present invention. The receiver has one or more receiving antenna, and the received signal 
is input to a multiple output equalizer 31 as a baseband digital received signal r m (k) (where m=1, 2. M) including 
one or more (M) signals. The received signal r m (k) is generated in a manner as shown in Fig. 30B, for example. 
[0141] The multiple output equalizer 31 is constructed in the same manner as shown in Fig. 2, and performs according 
to a processing procedure as shown in Fig. 3. Accordingly, an extrinsic information Mbfl)] is subtracted from a log- 
likelihood ratio A 2 [b(i)] from a decoder 24 shown in Fig. 22 in a subtracter 25, and the subtracted output is interleaved 
by an interleaver 26 to provide a priori information A^bG)], which is then converted in a series-parallel converter 15 
into N series of a priori information ^ [b^k)], [b N (k)] to be input to the multiple output equalizer 31 . 
[0142] Accordingly, N series of received signals are subject to a linear equalization in the multiple output equalizer 
31 in the similar manner as mentioned previously, delivering N log-likelihood ratio series A 1 [b^k)], — A 1 [b N (k)] t which 
are then input to a parallel-series converter 16 to be converted into a single log-likelihood ratio series A^bO)] to be 
supplied to a subtracter 22. With this arrangement, the input signal format to the multiple output equalizer 31 is similar 
to that described in connection with Figs. 1 to 3, and accordingly, N series log-likelihood ratios A^b^k)], - , [b N (k)J 
can be obtained by the equalization which has been mentioned above with reference to Figs. 1 to 3, and it will be 
readily seen that an iterative decoding processing is permitted by the use of the series-parallel converter 15 and the 
parallel-series conv rt r16. Inamann rcorr spending to a transmitted signal from a n-th transmitt r in the arrange- 
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ment of Figs. 1 to 3, an n-th or (n-th column) transmitted signal among N parall I transmitt d signals will be qualized. 
It will be readily seen that the mbbdiment described above in connection with Figs. 4 to 7 is applicable to th reception 
oftheparall I transmission of th Ns ri s signals. By processing in a plurality of qualizerstag swhichar incascad 
conn ction as illustrat d in Figs. 18 1 21. th reception characteristic can be improved as compared to the processing 
through a single equalizer stage as shown in Figs. 1 to 31. 

[0143] The turbo-reception method and the turbo-receiver according to the present invention are also applicable to 
the reception of convoluted code/turbo-code + interleaver + multi-value modulation such as QPSK 8PSK 16QAM 
64QAM etc.,, TCM (trellis coded modulation) / turbo TCM. 

Generation of M received signals 

[0144] M received signals ri (k), r M (k) are derived from M antennas #1 , ■» #M, but may be derived from a single 
antenna. Alternatively, M (which is greater than L) received signals may be obtained from L (which is an integer equal 
to or greater than 2) antennas. While not specifically shown in Fig. 1, received signals from antennas #1, #M are 
converted into baseband received signals ■« r M in a baseband converter and sampled to provide digital signals u 
(k). •". r M (k) at discrete time k. 1 
[0145] As shown in Fig. 30b, for example, received signals received by L=2 antennas #1 and #2 may be converted 
into baseband signals in baseband converters 61-1 and 61-2, respectively, the outputs of which may be sampled by 
A/D converters 64-1, 64-2 and 64-3, 64-4 in accordance with a sampling signal from a sampling signal generator 62 
and a sampling signal which is obtained by shifting the phase of the first mentioned sampling signal by one half period 
thereof T/2 in a phase shifter 63 to provide digital signals r^k), r 2 (k), r 3 (k), r 4 (k). respectively. Then the digital signal 
may be input to a turbo-receiver 30 as shown in Fig. 1 or Fig. 18 or Fig. 22 to provide N decoded outputs. It should be 
understood that the frequency of the sampling signal from the sampling signal generator 62 is chosen so that each 
sampling period of the received signals r^k), ■» r 4 (k) input to the turbo-receiver 30 coincides with the sampling period 
which is used when a single received signal r M (k) is received per antenna. 

EFFECTS OF THE INVENTION 



[0146] As discussed above, according to the first aspect of the invention (1 ), there is realized a multiple output (MIMO) 
reception method. To illustrate a quantitative effect, an error rate response is graphically shown in Figs. 23 and 24. In 
each Figure, E^ on the abscissa represents a bit power to noise ratio. A simulation presumes following conditions: 



Number of users (transmitters) 


2 


The number of multipaths of each user Q 


5 


Number of reception antennas 


2 


Number of information symbols per frame 


450 bits 


Number of unique words per frame 


25 bits 


Channel estimation method 


RLS (obliteration factor 0.99) 


Error correction code 


rate 1/2, constrained length 3 convoluted code 


Doppler frequency 


1000Hz (Rayleigh fading) 


Modulation 


BPSK 


Transmission rate 


20Mbps 


Decoder 24 


Max-Log-Map decoder 


Number of iterations 


4 


No fading within a frame 





No approximation by the matrix inversion lemma for the inverse matrix is used in the calculation of the filter coefficient w. 
[0147] Fig. 23 shows an error rate characteristic when the channel estimation is achieve perfectly (without no esti- 
mation error) or when the channels are known for a number of users N=2. a number of reception antennas M=2 a 
number of Raleigh paths Q=5. It will be seen that initial run is not an iteration and the second iteration is the first 
iteration. It will be seen that the error rate characteristic is significantly improved by the iteration. It will be seen from 
this that the turbo-reception method according to the invention for MIMO operates in a proper manner. 
[0148] Fig. 24 illustrates the effect of an iterative channel estimation (according to the fourth aspect of the invention) 
The abscissa indicates a threshold value Th. EJU Q is fixed to 4 dB (where E b is for one user), and Th=1.0 may be 
considered as a conventional method in which no symbol hard decision value is selected, or where a channel estimation 
using a symbol hard decision value is not made. In this instance, it will b s en from th drawing that BER characteristic 
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has a littl effect of iteration becaus the channel estimation is inaccurat . A threshold Th=0 indicates that all of hard 
decision valu s are directly used, and wh n hard d cision valu s of inf rmation symbols are used in this mann r, it 
willb apparent from the drawing that am an bit rrorrat isimprov d, allowing the channel stimationtobep rform' d 
xactly to a corresponding d gree. For a threshold value on the order of 0.2 to 0.6, the mean bit error rate is reduced 
than when Th=0 t indicating that it is preferable that only a hard decision value which is likely to be certain be used. It 
will be seen that Th which is around 0.25 is most preferable. 

[0149] Fig. 25 shows an error rate characteristic of an MIMO reception method in which a transmitted symbol hard 
decision value which is determined to be likely to be certain according to a threshold value is used in the channel 
estimation or which employs an iterative channel estimation in a form of a curve 66. In this instance, the threshold is 
set up at 0.25 and the result shown represents a characteristic after four iterations with parameters that N=2, M=2, 
Q-5Rayleigh, f d T s =1/20,000, and 900 symbols/frame. For purpose of comparison, an error rate characteristic with a 
perfect channel estimation is shown as a curve 67 while an error rate characteristic in which a hard decision value of 
information symbol is not used in a channel estimation or only a single channel estimation is made without iteration is 
shown by a curve 68. It will be seen from this graphic illustration that when the iterative channel estimation is used, 
the error rate characteristic approaches to that obtained by a perfect channel estimation. 

[0150] With the channel estimation method mentioned above, by determining whether a hard decision value is or is 
not likely to be certain on the basis of a decoded soft decision value, and by using symbol information having a hard 
decision value which is likely to be certain in the channel estimation during the next iteration, the channel estimation 
can be performed more correctly, allowing a decoding quality to be improved. 

[01 51 ] In order to confirm the effect of an embodiment in which a covariance matrix 0 (for noise other than Gaussian 
noise) is estimated, a simulation is made with parameters indicated below. 



Number of users (transmitters) 


3 (one of which is chosen as an unknown interference: i(k)) 


The number of multipaths of each user Q 


5 


Number of reception antennas 


3 


Number of information symbols per frame 


450 bits 


Error correction code 


rate 1/2, constrained length 3 convoluted code 


Doppler frequency 


1000Hz (Rayleigh fading) 


Modulation 


BPSK 


Transmission rate 


20Mbps 


Decoder 24 


Log-MAP decoder 


Number of iterations 


4 



[01 52] Three users (transmitters) are chosen to be of an equal power. 

[0153] Fig. 26 shows a result of simulation of BER (bit error rate) characteristic of a turbo-receiver which estimates 
H and 0 as shown in Figs. 14, 15 and 16, and Fig. 27 shows BER characteristic which directly uses the turbo-receiver 
shown in Fig. 1 (which employs the method shown in Fig. 13). In Fig. 26, it is assumed that the noise comprises only 
white Gaussian noise, and it is seen that two or more iterations of the channel estimation and the decoding processing 
has little effect. By contrast, it will be seen from Fig. 27 that as the number of iterations is increased, an improvement 
in the BER characteristic is achieved, and in addition, BER is considerably less than that shown in Fig. 26 for the same 

[0154] In order to confirm the effect of the embodiment (according to the second aspect of the invention) in which 
an error correction decoding result is reflected into a symbol soft decision value b' n (k) of a received signal from an 
intended user (transmitter), a simulation is made with parameters as indicated below 



Number of users (transmitters) N 


4 


The number of multipaths of each user Q 


5 


Number of reception antennas 


2 


Number of information symbols per frame 


900 bits 


Error correction code 


convoluted code (coded rate:1/2, constrained length 3) 


Modulation 


BPSK 


Decoder 24 


Log-Map decoder 


Erroneous coded rate 


1/2 


Number of iterations 


5 
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10 



15 



20 



25 



30 



[01 55] A choice is made that f(b' n (k))= a x b' n (k) 

[0156] Fig. 28 shows BER characteristic of a multiple output turbo-receiver shown in Fig. 1 and a multiple input 
multipl output turbo-receiver in which an rror correction decoding r suit is reflect d into b' n (k), with plott d points 
shown in black for th former and in whit for th latt r. It is to be noted that in the plot, a circle represents an initial 
run, a down-directed triangle a second iteration, a diamond a third iteration, a left-directed triangle a fourth iteration 
and a right-directed triangle a fifth iteration. Fig. 28A shows a result of a simulation for the BER characteristic plotted 
against E^ when a is fixed to 0.2, and Fig. 28B shows a result of simulation for the BER characteristic plotted against 
a when Eb/N 0 is 6dB. It will be seen that a=0 represents that b' n (k)=0. It is noted from Fig. 28A that in the multiple input 
multiple output receiver in which an error correction decoding result is reflected into b* n (k). an improvement in BER 
after a third and a subsequent iteration is greater in comparison to BER obtained during a preceding iteration, as 
compared with the multiple input multiple output turbo-receiver shown in Fig. 1 , and that after a third and a subsequent 
iteration, when EJN Q required to achieve each BER in a range BER>10r 4 is compared, the multiple input multiple 
output turbo-receiver in which an error correction decoding result is reflected into b' n (k) has a gain of about 0.5 dB or 
greater in comparison to the multiple input multiple output turbo-receiver shown in Fig. 1. It is also noted that during 
the fifth iteration at Eb/N 0 =6dB, BER=10-5 is attained, which is a reduction in BER by a factor of 1/10 or more in 
comparison to the receiver shown in Fig. 1. It will be seen from Fig. 28B that an improvement is achieved in a range 
of a indicated by an equality 0«x<0.6 and that when a exceeds 0.6, the BER characteristic becomes degraded, pre- 
venting a correct decoding result from being obtained. From this result, it is seen that an optimum value of a is 0.2. 
However, it should be understood that the value of a is not limited to the optimum value, but a suitable range of a 
having an improving effect can be changed depending on the number of users to be received, a propagation environ- 
ment having interferences, a number of receiving antennas or the like. In addition, a different value may be chosen as 
an optimum value for a. 

[0157] When BPSK modulation is usedfora number of users (transmitters) which is equal to N, a number of multipaths 
from each transmitter which is equal to Q and a number of receiver antennas which is equal to M, the quantity of 
calculation which is required in an equalizer when a conventional single user turbo-receiver is directly extended to a 
multiple output (MIMO) is on the order of 2 N (Q- 1 ), as mentioned previously, but with the turbo-reception method accord- 
ing to the third aspect of the invention, the quantity of calculation can be reduced to the oreder of N(MQ) 3 By way of 
example, assuming that N=8, Q=20and M=8, 2N(Q-D«5-1 0*5 while N(MQ)3«37-107 thus demonstrating that the quantity 
of calculation can drastically be reduced according to the turbo-reception method according to the second aspect of 
the present invention. 

[0158] A simulation has been conducted under conditions given below in order to confirm that a good bit error rate 
characteristic can be obtained according to the turbo-reception method according to the third aspect of the invention. 
It is assumed that the channel matrix H is known. 
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Number of users N 


4 


The number of multipaths of each user Q 


5 


Number of reception antennas 


2 


Number of information symbols per frame 


900 bits 


Error correction code 


rate: 1/2, constrained length 3 convoluted code 


Doppler Frequency 


1000Hz (Rayleigh fading) 


Modulation 


BPSK 


Decoder 24 


Log-Map decoder 


Transmission rate 


20Mbps 


Decoder 


Log-MAP decoder 


Number of iterations 


6 



[01 59] It is assumed that the channel estimation takes place in an ideal manner. 

[01 60] Fig. 29 shows a result of simulation for BER (bit error rate characteristic). The abscissa represents mean Eb 
(bit power)/ N 0 (noise power), denotation f d in the graph represents a doppler frequency and D s the period of a trans- 
mitted symbol. MRC shown on this graph is BER characteristic upon Viterbi decoding of a signal after a maximal ratio 
combining on an order 10 (two antennasxS paths) diversity channels, and this corresponds to the BER characteristic 
which results when the equalizer has perfectly canceled interferences. Thus, the quality of the receiver can be evaluated 
by seeing how the BER after the iterations is located close to the MRC curve. It is seen from Fig. 29 that with the turbo- 
reception method according to the second aspect of the present invention, the higher E^Nq, the less BER, and the 
greater the number of iterations, the more the BER characteristic approaches MRC BRE characteristic, in particular 
BER characteristic is located very close to MRC with a number of iterations qual to 6. Thus, it is confirmed that th 
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muttipl output turbo-receiv r which employs the turbo-reception method according to the third asp ct f the present 
invention op rat s properly under severe conditions of four users, each having fiv paths and using two reception 
antennas. 



5 

Claims 



10 



15 



1. A turbo-reception method of receiving a signal from N transmitters where N is an integer equal to or greater than 
2, comprising the steps of 

calculating a channel value h mn (q) and a channel matrix H from M received signals r m (where M is an integer 
equal to or greater than 1) and a known signal, where m=1, M; n=1, N and q=0, Q-1, and Q represents 
a number of multipaths of each transmitted wave; 

determining a soft decision transmitted symbol b' n (k) from N a priori information ^ [b n (k)] where k represents 
a discrete time; 

using the channel value h mn (q) and the soft decision transmitted symbol b' n (k), calculating an interference 
component H-B'fk) upon a transmitted signal from an n-th transmitter, 
where 



20 



25 



H( 0 ) — H (Q _j) 0 
0 H( 0 ) — Hflj.,) 



30 



H (q) = 



MKq) 



h lN(q) 



35 



40 



B , (k)=[b fT (k+Q-1)...b' T (k)...b ,T (k-Q+1)] T 
b'(k+q)=[b f 1 (k+q)b' 2 (k+q)...b' N (k+q)] T 
q=Q-1 -Q+1 q*0 



b'(k)=[b' 1 (k)...O...b' N (k)] T q=0 

b'(k) having a zero element at an n-th position and [ ] T representing a transposed matrix; 

subtracting the interferenced component H-B'(k) from a received matrix y(k) to obtain a difference matrix y' 
(k) where 



50 y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 

r(k)=[ ri (k)r 2 (k)...r M (k)] T 

55 usin 9 foe channel matrix H or a reference signal to determine an adaptive filter coefficient w n (k) which is to 

be applied to a received signal corresponding to a transmitted signal from an n-th transmitter in order to eliminate 
residual interference components in the difference matrix y'(k); 

andfilt ring the diff rence matrix y'(k) with th adaptiv filter co fficient w n (k) to provide a log-lik lihood ratio 
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for th received signal which corr sp nds to the transmitted signal from th n-th transmitter and from which int r- 
f rences are eliminated. 

2. A turbo-reception method according to Claim 1 in which d noting a covariance matrix of noise components in the 
received matrix y (k) by U, using the soft decision transmitted symbol b' n (k) and the channel matrix H to calculate 
the adaptive filter coefficient w n (k) according to the following formula: 

w n (k)=(HG(k)H H +U)' 1 h 
G(k)=diagtD(k+Q-1)...D(k)...D(k-Q+1)] 



15 



20 



DOc+q^agfl-b'^Ck+q),. . .,1-b' 2 n (k+q),. ..,1-b' 2 N(k+q)] 
q=Q-l...-Q+l,q#) 
=diag[l-b ,2 ,(k+q),...,l,...,l-b ,2 N (k+q)] q=0 



25 



30 



h = 



H l.(Q-l).tf+n 
H 2.(Q-l).N+n 



H, 



M.(Q-l).H+n 



35 



40 



45 



where H^o.^.^ represents an element of the matrix H located at row 1 and column (Q-1)N+n. 

A turbo-reception method of receiving a signal from N transmitters where N is an integer equal to or greater than 
2, comprising the steps of 

calculating a channel value h mn (q) and a channel matrix H from M received signals r m where M is an integer 
equal to or greater than 1 and where m is 1, — ,M; n=1 , N; q=0, ■», Q-1, and Q represents a number of multipaths 
of each transmitted wave; 

determining a soft decision transmitted symbol b' n (k) from N a priori information ^ [b n (k)] where k represents 
a discrete time; 

using the channel value h mn (q) and the soft decision transmitted symbol b' n (k) to calculate an interference 
component H-B'(k) upon a transmitted signal from the n-th transmitter where 



50 



H= 



H( 0 ) - tyg.j) 0 

o H( 0 ) ••- H(Q_j) 
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ll(q) — h lN(q> 



h Ml(q) htfN(q) 



B'(k)=[b ,T (k+Q-1)...b' T (k)...b' T (k-Q+1)] T 



b , (k+q)=[b' 1 (k+q)b f 2 (k+q)...b' N (k+q)] T 
q=Q-1 -Q+1 q*0 



b , (k)=[b' 1 (k)...-f(b' n (k))...b' N (k)] T q=0 

b'(k) having an element f(b' n (k)) at an n-th position, f ( ) being a function of a variable b' n (k) and satisfying that f(0) 
= 0 and 

d{f(b' n (k))}/d{b , n (k)}^0, and [ ] T representing a transposed matrix; 

subtracting the interference component H-B'(k) from a received matrix y(k) to determine a difference matrix 
y'(k) where 



y(k)=[r T (k+Q-1)r T (k+Q-2)... r T (k)] T 



r(k)=[ ri (k)r 2 (k)...r M (k)] T 

using the channel matrix H or reference signal to determine an adaptive filter coefficient w n (k) which is to be 
applied to a received signal corresponding to the transmitted signal from the n-th transmitter in order to eliminate 
residual interference components which remain in the difference matrix y'(k) ; 

and filtering the difference matrix y'(k) with the adaptive filter coefficient w n (k) to provide a log-likelihood ratio 
for the received signal which corresponds to the transmitted signal from the n-th transmitter and from which inter- 
ferences are eliminated. 

A turbo-reception method according to Claim 3, further comprising, denoting a covariance matrix of noise compo- 
nents in the received matrix y(k) by U, using the soft decision transmitted symbol b' n (k) and a channel matrix H to 
calculate the adaptive filter coefficient w n (k) according to the following formula: 

w n (k)=(HG(k)H H +U)* 1 h 



G(k)=diag[D(k+Q-1)...D(k)...D(k-Q+1)] 



H (q) = 
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D(k+q)^ag[l-b^(k+q),...a-b'\(k+q),...4-b ,2 N (k4< 1 )] 
q=Q-l...-QH,q#) 

^ag[l-b*(k+q),...,l-b' 2 „,^ 
l-b ,2 n+I (k),...,l-b ,2 N (k+q)] q=0 



where E[ ] represents a mean value, and 



^1 . (Q-1) . N+n 
H 2.(Q-l).N+n 

L H M,(Q-l).N+n_ 



where H 1 ^Q_ 1 j. N+n represents an element of the matrix H at row 1 and a column (Q-1) N+n. 

5. A turbo-reception method according to Claim 2 or 4 in which an inverse matrix calculation during the calculation 
of the adaptive filter coefficient w n (k) takes place by using matrix inversion lemma. 

6. A turbo-reception method according to Claim 1 or 2 in which the covariance matrix U of noise components in the 
received matrix y (k) is defined as a 2 I which is determined by a variance a 2 of a Gaussian distribution and a unit 
matrix. 

7. A turbo-reception method according to Claim 1 or 2 in which the covariance matrix U of noise components in the 
received matrix y (k) is derived, using the received matrix y(k) and the estimated channel matrix H, as follows: 

0= 2k=o Tr (y(k)-H B(k)) (y(k)-H B(k)) H 
B(k)=[b T (k+Q-1)...b T (k)...b T (k-Q+1)] T 



b(k+q)=[b 1 (k+q)...b N (k+q)] T 
(q=-Q+1...Q-1) 

where Tr represents the length of the reference signal. 

8. A turbo-reception method according Claim 2 or 3 in which D(k+q) is approximated by 0 for q*0, and D(k+q) is 
approximated by diag [0, •-, 1 , »-, 0] for q=0. 

9. A turbo-reception method of receiving a signal from N transmitters where N is an integer equal to or greater than 
2, comprising the steps of 

calculating a channel value h mn (q) and a channel matrix H from M received signals r m (where M is an integer 
qual to or greater than 1) and a known signal wh re m=1, ••, M; n=1, -, N; q=0, Q-1, and Q repres nts a 
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number of multipaths f each transmitted wave; 

d t rrniningasoftd cision transmitt d symbol b* n (k) from N a priori information X2[b n (k)] where k rep res nts 
a discrete time; 

using th chann I value h mn (q) and the soft decision transmitted symbol b' n (k) to calculate an interference 
component H-B f (k) upon a transmitted signal from an n-th transmitter where 



H= 



H( 0 ) ••• Hflj.j) o 

* * • 

0 tyo) — H (Q-1) 



H (q) = 



l l 1 (q) 



h lN(q) 



B f (k)=[b ,T (k+Q-1)...b ,T (k)...b ,T (k-Q+1)] T 

b , (k+q)=[b I 1 (k+q)b' 2 (k+q)...b' N (k+q)] T 
q=Q-1 -Q+1 q*0 



b , (k)=[b , 1 (k)...O...b' N (k)] T q=0 

b'(k) having a zero element at an n-th position and Q T representing a transposed matrix; 

subtracting the interference component H-B'(k) from a received matrix y(k) to determine a difference matrix 
y'(k) where 

y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 
r(k)=[r 1 (k)r 2 (k)...r M (k)] T 

defining a covariance matrix of noise components in the received matrix y(k) as a 2 I which is determined 
from a variance a 2 of a Gaussian distribution and a unit matrix I, and filtering the difference matrix y'(k) with an 
adaptive filter coefficient w n which is determined by 
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10 



h = 



H l.(Q-l).N+n 
^2.(Q-l).N+n 



H 



M-Q.(Q-l).N+n 



15 



20 



25 



to provide a log-likelihood ratio for the received signal which corresponds to the transmitted signal from the n-th 
transmitter and from which interferences are eliminated. 

10. A turbo-reception method of receiving a signal from N transmitters where N is an integer equal to or greater than 
2, comprising the steps of 

calculating a channel value h nm (q) and a channel matrix H from M received signals r m (where M is an integer 
equal to or greater than 1) and a known signal where m=1, ■», M ; n=1, N; q=0, Q-1, and Q represents a 
number of multi paths of each transmitted wave; 

determining a soft decision transmitted symbol b* n (k) from N a priori information [b n (k)] where k represents 
a discrete time; 

calculating an interference component H-B'(k) upon a transmitted signal from an n-th transmitter using the 
channel value h mn (q) and the soft decision transmitted symbol b' n (k) as follows: 



30 



H = 



H(0) 
o 



H (Q-1) 

H (o) 



o 



35 



h ll(q) — h lN(q) 



h Ml(q) — tMN(q) 



B , (k)=[b ,T (k+Q-1)...b' T (k)...b ,T (k-Q+1)] T 

45 

b'(k+q)=[b , 1 (k+q)b' 2 (k+q)...b , N (k+q)] T 
q=Q-1 -Q+1 q*0 

so 

b'(k)=[b' 1 (k)...-f(b' n (k))...b' N (k)] T q=0 



where b'(k) has an element f(b' n (k)) at an n-th position, f( ) is a function of a variable b' n (k) which satisfies that of 
f(0) =0 and is d{f(b' n (k))}/d{ b' n (k)} a 0 and [ ] T is a transposed matrix; 

subtracting the interference component H-B'(k) from a received matrix y(k) to define a difference matrix y* 
(k) where 



H (q) = 

40 
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y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 



r(k)=[ ri (k)r 2 (k)...r M (k)] T 

defining a covariance matrix for noise components in the received matrix y(k) as a 2 I which is determined 
from a variance a 2 of a Gaussian distribution and unit matrix I, and filtering the difference matrix y'(k) with an 
adaptive filter coefficient w n (k) which is determined by 



H l.(Q-l).N+n 
H 2.(Q-l).N+n 



u H M.Q.(Q-l).N+n 



to provide a log-likelihood ratio for the received signal which corresponds to the transmitted signal from the n-th 
transmitter and from which interferences are eliminated. 

11. A turbo-reception method of receiving a transmitted signal from N transmitters where N is an integer equal to or 
greater than 2 comprising the steps of 

determining a channel value which is a transmission characteristic of a received signal from M received 
signals (where M is an integer equal to or greater than 1) and a known signal;. 

estimating a soft decision transmitted symbol each from N a priori information; 

dividing N transmitted signals each into L (L£N) groups of transmitted signals including one or a plurality of 
transmitted signals, using a soft decision transmitted symbol and a channel matrix comprising channel values to 
determine L equalized signals from which interferences from other group of transmitting signals are eliminated 
and post-equalization channel information each corresponding to the transmission characteristic of the equalized 
signal for each group of transmitted signals; 

for each combination of L equalized signals and associated channel information, treating the group of equal- 
ized signals as a received signal having a channel value defined by the channel information and where there is a 
plurality of transmitted signals which constitute the group of equalized signals, dividing such transmitted signal 
into a plurality of sub-group of one or a plurality of transmitted signals, using the soft decision transmitted signal 
to determine an equalized signal for the subgroup of transmitted signals from which interferences from the other 
groups of transmitted signals are eliminated and associated post-equalization channel information, and where the 
group comprises a single transmitted signal, using the resulting equalized signal, the channel information and the 
soft decision transmitted signal to eliminate an interference of the transmitted signal itself due to multipaths; 

repeating the steps of dividing, eliminating interferences and generating post-equalization channel informa- 
tion until transmitted signals which constitute every equalized signal becomes a single one, thus finally determining 
a equalized signal for each transmitted signal from which an interference caused by multipaths of its own is elim- 
inated or determining, for the combination of the equalized signal and associated channel information, an equalized 
signal from which an interference between different transmitted signals and an intersymbol interference by the 
transmitted signal itself is eliminated for each transmitted signal which constitutes the equalized signal. 

12. A turbo-reception method according to Claim 11, in which for each group of transmitted signals, using the soft 
decision transmitted symbol and the channel value to generate a replica of interferences caused by other groups 
of transmitted signals, subtracting the replica of interferences from the received signal to define a difference signal, 
determining, for each difference signal, a filter characteristic which is used to eliminate residual interference and 
associated post-equalization channel information from the channel value and the soft decision transmitted symbol, 
filtering a corresponding difference signal with the residual interference eliminating filter characteristic to derive 
th qualiz d signal. 



h = 
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13. A turbo-reception method according to Claim 12, comprising the steps of 

determining from th receiv d signals r^k), r M (k), a receiv d matrix as defined below 



y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 



r(k)=[ ri (k)r 2 (k)...r M (k)] T 

where [ ] T represents a transposed matrix; 

defining the transmission characteristic as a channel matrix H as defined below 



H= 



H(o> — H (Q _ 0 o 

• ■ 

o H( 0 ) — H(q-d 



H(q)=[h 1 (q)...h N (q)] 
h n (q)=[h 1n (q)...h Mn (q)] T 

where m=1 t M; n=1, — , N; and q=0, — , Q-1 and Q is a number of mutipaths, and h mn (q) is a channel value of 
a path q from an n-th transmitter which is contained in the received signal r m ; 

denoting the soft decision transmitted symbol by b' n (k), calculating a replica H-B'(k) upon one of groups of 
transmitted signals, formed by a first to a U-th transmitted signal where U is an integer satisfying N>U2s1 caused 
by other groups of transmitted signals as follows; 

B'(k)=[b ,T (k+Q-1)...b' T (k)...b' T (k-Q+1)] T 

b'tk+qMb*., (k+q)b' 2 (k+q)...b' n (k+q)...b' N (k+q)] T 
:q=Q-1,...1 

b'(k+q)=[0...0 b u+1 '(k+q)...b' N (k+q)] T :q=0 -Q+1 

where b'(k+q) have zero elements which are equal to U in number; 

subtracting the interference replica H-B'(k) from the received matrix y(k) to define the difference matrix y* g (k). 

14. A turbo-reception method according to Claim 11 in which when a further elimination of interferences is desired for 
the equalized signal and associated channnel information, the number of multipaths which are used during the 
elimination of interferences from the equalized signal is reduced. 

15. A turbo-reception method according to Claim 14, further comprising the steps of determining from the received 
signals ^(k), r M (k), a received matrix defined as follows; 



y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 
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r(k)=[r 1 (k)r 2 (k)...r M (k)] T 

wher [ ] T r pres nts a transposed matrix; 

defining the transmission characteristic in the form of a channel matrix H, defined as follows; 



H= 



H(o) — H ( q_d o 
o H( 0) Hflj.D 



H(q)=[h 1 (q)...h N (q)] 



h n (q)=[h 1n (q)...h Mn (q)] T 

where m=1, -, M; n=1 t — , N; and q=0, — , Q-1 and Q is a number of muti paths, and h mn (q) is a channel value of 
a path q from an n-th transmitter which is contained in the received signal r m ; 

denoting the soft decision transmitted symbol by b* n (k), one of the groups of transmitted signals comprising 
a first to U-th transmitted signal where U is an integer satisfying the inequality N>lte1, an equalized signal being 
formed for this group of transmitted signals so that interferences therein are eliminated by considering a number 
of multipaths equal to Q* where Q'<Q, calculating a replica of interference upon this group from other groups of 
transmitted in the form of H-B'(k) as follows; 



B'(k)=[b ,T (k+Q-1)...b ,T (k)...b' T (k-Q+1)] T 

b , (k+q)=[b , 1 (k+q)b' 2 (k+q)...b , n (k+q)...b t N (k+q)] T :q=Q-1 ,...1 

b'(k+q)=[0...0 b' u+1 (k+q)...b' N (k+q)] T :q=0 -Q'+1 

where b'(k+q) have 0 elements which are equal to U in number, and 

b , (k+q)=[b , 1 (k+q)...b' n (k+q)...b' N (k+q)] T :q=Q' -Q+1 

subtracting the interference replica H-B'(k) from the received matrix y(k) to define the difference matrix y' g (k). 

16. A time-reception method according to one of Claims 1, 3, 9, 10 and 11 in which during a second and a subsequent 
iteration of the turbo-reception processing, both the known singnal and the transmittes symbol hard decision output 
which is obtained during a previous alteration are used as reference signals, and the reference signals and the 
received signal are used to calculate the channel matrix. 

17. A turbo-reception method according to Claim 16 in which one of the transmitted symbol hard decision outputs 
which are obtained during a previous iteration and which has a certainty in excess of a given value is also used 
as a reference signal to be used in the calculation of the channel matrix. 

18. A turbo-reception method according to one of Claims 1, 3, 9, 10 and 11 in which N a priori information ^ [b n (k)] 
are derived from N decoders which correspond to the N transmitters, a log-likelihood ratio which is obtained for 
the received signal which corresponds to the n-th transmitted signal and through which interferences are eliminated 
is fed to a corresponding one of the decoders. 
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9. A turbo-reception m thod according to on of Claims 1 t 3, 9, 10 and 11 in which th N transmitt d signals are 
transmitted from N transmitters which transmit a single information sen s in the form of N parallel s ries, the N a 
priori informati n information [b n (k)] being a result of sen s-parallel converstion of a priori information [bO)l 
from one of th d coders, N log-lik lihood ratios which are received signals corresponding t th N transmitted 
signals and from which interferences are eliminated being subject to a parallel-series convertion before being fed 
to the decoders. 

0. A turbo-receiver for receiving a signal from N transmitters where N represents an integer equal to or greater than 
2, comprising 

a received signal generator for forming M received signals r m (where M is an integer equal to or greater than 
1) where m=1, M; 

a channel estimator to which each received signal r m and a reference signal in the form of a known signal 
are input for calculating a channel value h^q) and a channel matrix H, as defined below 



H = 



H(o> — H (Q _ 0 o 
0 H (0 ) ... H (Q _ 0 



h ll(q) — h lN(q) 

• • ♦ 

• • • 

• • ■ 

h Ml(q) — hMN(q)_ 

n=l, ... N 

a received matrix generator which receives each received signal r m for generating a received matrix, as 
defined below 

y(k)=[r T (k+Q-1) r T (k+Q-2)... r T (k)] T 



r(k)=[ ri (k)r 2 (k)...r M (k)] T 

where k represents a discrete time, Q a number of multipaths of each transmitted wave, q=0, Q-1, and [ ] T 
representing a transposed matrix; 

a soft decision symbol generator receiving N a priori information for generating a soft decision transmitted 
symbol b' n (k); 

a replica matrix generator to which respective soft decision transmitted symbols b^k) to b N '(k) are input to 
generate an interference replica matrix B'(k) with respect to an n-th transmitted signal as follows; 

B'(k)=[b ,T (k+Q-1)...b ,T (k)...b ,T (k-Q+1)] T 



b'(k+q)=[b^ (k+q)b , 2 (k+q)...b' N (k+q)] T 
q=Q-1 -Q+1, q*0 



H (q) = 
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b'fk)^ (k)...O...b' N (k)] T q=0 

where b'(k) has a z ro element at n-th position; 

a filter processor to which the channel matrix H and the interference replica matrix B'(k) are input for calcu- 
lating and delivering an interference component H-B'(k) with respect to a received signal corresponding to the n- 
th transmitted signal; 

a difference calculator to which the interference component H-B'(k) and the received matrix y(k) are input 
for delivering a difference matrix y'(k) = y(k) - H-B'(k); 

a filter coefficient estimator to which the channel matrix H or a reference signal is input for determining an 
adaptive filter coefficient w n (k)to be applied to a received signal which corresponds to a transmitted signal from 
the n-th transmitter in order to eliminate residual interference components which remain in the reference matrix y' 
(k); 

and an adaptive filter to which the difference matrix y'(k) and the adaptive filter coefficient w n (k) are input to 
filter y'(k) to provide a log-likelihood ratio as an received signal which corresponds to the transmitted signal from 
the n-th transmitter and from which interference is eliminated and which is then fed to the n-th decoder. 

21 . A turbo-receiver for receiving a signal from N transmitters where N is an integer equal to or greater than 2, com- 
prising 

a received signal generator for forming M received signals r m (where M is an integer equal to or greater than 
1) where m=1, — t M; 
N decoders; 

a channel estimator to which each received signal r m and a reference signal in the form of a known signal 
are input to calculate a channel value h mn (q) and a channel matrix H, as defined below. 



H( 0 ) — H (Q _ 0 o 
0 H( 0 ) — H( Q _ j) 



h ll(q) — h lN(q) 

■ • • 

• • • 

• * • 

h Hl(q) — hfflKq). 

n= 1, ...,N 

a received matrix generator to which each received signal r m is input to generate a received matrix, as defined 
below 

y(k)=[r T (k+Q-1)r T (k+Q-2)...r T (k)] T 



r(k)=[r 1 (k)r 2 (k)...r M (k)] T 

where k represents a discrete time. Q a number of multipaths of each transmitted wave, q=0. —., Q-1 and [ ] T a 
transposed matrix; 

a soft decision symbol generator to which N a priori information are input to generate soft decision transmitted 
symbols b' n (k) (where n=1, N); 

a replica matrix generator to which th soft decision transmitted symbols b^ (k) to b' n (k) are input to gen rate 
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an int rf rence replica matrix B'(k) for a transmitted signal from an n-th transmitt r as indicat d below 

B'(k)=[b' T (k+Q-1)...b ,T (k)...b ,T (k-Q+1)] T 



b , (k+q)=[b , 1 (k+q)b' 2 (k+q)...b' N (k+q)] T 
q=Q-1 -Q+1 t q^0 

b t (k)=[b' 1 (k)...-f(b' n (k))...b' N (k)] T q=0 

where b'(k) has an element f(b* n (k)) at an n-th position, and f ( ) is a function of a variable b' n (k) which satisfies 
that f(0)=0 and 



d{f(b' n (k))}/{b' n (k)}2>0; 

a filter processor to which the channel matrix H and the interference replica matrix B'(k) are input to calculate 
and deliver an interference component H-B'(k) for a received signal which corresponds to the transmitted signal 
from the n-th transmitter; 

a difference calculator to which the interference component H-B*(k) and the received matrix y(k) are input to 
deliver a difference matrix y'(k) =y(k)- H-B'(k); 

a filter coefficient estimator to which the channel matrix H or a reference signal is input to determine an 
adaptive filter coefficient w n (k) for a received signal which corresponds to a transmitted signal from the n-th trans- 
mitter in order to eliminate residual interferences in the difference matrix y'(k); 

and an adaptive filter to which the difference matrix y'(k) and the adaptive filter coefficient w n (k) are input to 
filter y'(k) to provide a log-likelihood ratio as a received signal which corresponds to the transmitted signal from 
the n-th transmitter and from which interferences are eliminated and which is then fed to the n-th decoder. 

22. A turbo-receiver for receiving a transmitted signal from N transmitters where N is an integer equal to or greater 
than 2, comprising 

a received signal generator for forming M received signals (where M is an integer equal to or greater than 1 ); 

a channel estimator to which the N received signals and a reference signals in the form of a known signal 
are input to estimate a channel value representing a transmission characteristic; 

a prestage equalizer to which the M received signals, the channel value and N a priori information are input 
for delivering, for each transmitted signal from one or more transmitters, a plurality of sets of equalized signal from 
which interferences by transmitted signals from other transmitters are eliminated and an associated post-equali- 
zation channel information; 

and a plurality of poststage equalizers to which the set of the equalized signal and associated channel in- 
formation from the prestage equalizer, and a priori information which corresponds to the transmitted signal which 
constitutes the equalized signal are input to deliver a log-likelihood ratio by eliminating from the equalized signal 
an intersymbol interference from each transmitted signal which constitutes the equalized signal due to the multiple 
paths and a mutual interferences between each transmitted signal and other transmitted signals which also con- 
stituted the equalized signal. 

23. A turbo-receiver for receiving a transmitted signal from N transmitters where N is an integer equal to or greater 
than 2 , comprising 

a received signal generator for forming M received signals (where M is an integer equal to or greater than 1 ); 

a channel estimator to which the N received signals and a reference signal in the form of a known signal are 
input to estimate a channel value representing a transmission characteristic; 

a prestage equalizer to which the M received signals, the channel value and N a priori information are input 
to deliver, for each transmitted signal from one or more transmitters, a plurality of sets of an equalized signal from 
which interferences by transmitted signals from other transmitters are eliminated and an associated post-equali- 
zation channel information; 

and a plurality of poststage equalizers to which the sets of the equalized signal and associated channel 
information from the prestage equalizer, and a priori information which corr sponds to a plurality of transmitted 
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signals which constitute the qualiz d signal are input to deliv r t for each r for a sub-group of transmitted signals 
among the plurality of transmitted signals which constituted the equalized signals, a plurality of sets of an qualized 
signal from which interference by oth r transmitt d signals which constitut d that the qualized signals are lim- 
inated and an associated post-equalization channel information. 

24. A turbo-receiver according to one of Claims 20 to 23, further comprising a previous symbol memory in which a 
hard decision transmitted symbol from the decoder is stored to be updated, and means for reading the hard decision 
transmitted symbol from the previous symbol memory to feed it to the channel estimator as the reference signal 
during a second and a subsequent iteration of the turbo-reception processing. 

25. A turbo-receiver according to Claim 24 further comprising a comparator for comparing a soft decision transmitted 
symbol which is input thereto against a threshold value, and a selector which is controlled by an output from the 
comparator so that one of the hard decision transmitted symbols, for which a corresponding soft decision trans- 
mitted symbol has a value in excess of the threshold value, is stored in the previous symbol memory. 

26. A turbo-receiver according to one of Claims 20 to 23, further comprising N decoders, to which the N log-likelihood 
ratios are delivered, the N decoders providing the N a priori information at their outputs. 

27. A turbo-receiver according to one of Claims 20 to 23 in which the N transmitted signals are transmitted signals as 
the N transmitters transmit a single information series as N parallel series, further comprising a parallel-series 
converter for converting the delivered N log-likelihood ratios into serial series, a decoder to which the log-likelihood 
ratio in the serial series is input, and a series-parallel converter for converting a priori information from the decoder 
into N parallel series to provide the N a priori information. 

28. A turbo-reception method in which a channel value representing a transmission path characteristic of a received 
signal is estimated from the received signal and a known signal serving as a reference signal, the received signal 
is processed in accordance with the estimated channel value, the processed signal is decoded, and the processing 
which uses the estimated channel value and the decoding are iterated upon the same received signal; further 
comprising 

determining the certainty of a decoded hard decision information signal on the basis of a value of a corre- 
sponding soft decision information symbol, and using a hard decision information symbol having a certainty which 
is in excess of a given value in the channel estimation of a next iteration as a reference signal. 

29. A reception method according to Claim 28, further comprising the step of calculating a 2 I as a covariance matrix 
of noise component in a received matrix y(k) where a 2 represents a variance of Gaussian distribution and I a unit 
matrix. 

30. A reception method according to Claim 28 further comprising using an estimated channel matrix H and a received 
signal matrix y(k) to calculate a covariance matrix U of noise components within the received signal matrix y(k) 
during each iteration as follows; 

0=Ik=0 Tr (y«-H'B(k))-(y(k)-H-B(k)) H 



B(k)=[b T (k+Q-1)...b T (k)...b T (k-Q+1)] T 

b(k+q)=[b 1 (k+q)...b N (k+q)] T (q=-Q+1...Q-1) 

bi (k+q) to b N (k+q) being the reference signals comprising the known signal and the hard decision information 
symbols having a certainty which is in excess of a given value, Tr being the length of the reference signal. 

31 . A reception method according to Claim 28 in which the iteration of the processing which uses the estimated channel 
value and decoding processing comprises an iteration of determining a linear equalization filter in accordance with 
the estimated channel value, processing the received signal by the linear equalization filter and decoding the 
processed signal. 
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32. A reception method according to Claim 28 in which the it ration of the processing which us sthe stimated channel 
value and decoding processing compris s an it ration of performing a rak synthesis which comp nsat s for a 
phas rotation which each symbol has expert need on a transmission path in a rak synthesis processor in ac- 
cordance with the estimat d channel valu , and decoding th signal which is produced by the rake synthesis in a 
turbo-decoder. 

33. A reception method according to Claim 28 in which the iteration of the processing which uses the estimated channel 
value and decoding processing comprises an iteration of setting up weights which determine an antenna directivity 
response of an adaptive array antenna receiver in accordance with the estimated channel value, and decoding an 
output from the adaptive array antenna receiver in a turbo-decoder. 

34. A reception method according to Claim 33 in which a rake synthesis which compensates for a phase rotation which 
each symbol has experienced on a transmission path in accordance with the estimated channel value is made for 
an output from the adaptive array antenna receiver in a rake synthesis processor, and a resulting signal from the 
rake synthesis is fed to the turbo-decoder. 

35. A receiver in which a channel value representing a transmission path characteristic of a received signal is estimated 
from the received signal and a known signal serving as a reference signal, the received signal is processed by 
using the estimated channel value, the processed signal is decoded, and the processing which uses the estimated 
channel value and the decoded processing are iterated upon the same received signal; further comprising 

means for determining whether or not the certainty of a decoded hard decision information symbol is in the 
excess of a given value by seeing whether or not a corresponding soft decision information symbol has a value 
which is in excess of a threshold value; 

and a previous symbol memory having its stored content updated by a hard decision information symbol 
which has been determined as likely to be certain, the stored content of the previous symbol memory being used 
as a reference signal during the channel estimation of a next ilteration. 
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